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ABSTRACT 


An experimental study was carried out to obtain the 
moisture migration characteristics of Devon Silt under 
varying conditions of applied load and freezing rate. All 
testing was carried out under open system drainage 


conditions. 


Short term and long term experiments were carried out 
and it was found that for the majority of cases short tern 
testing could be used to predict long term freezing 
behaviour. The limitations of this observation are not 
completely understood however guidelines as to duration of 
testing to obtain the true moisture migration 


characteristics are given. 


The experimental program showed that the moisture 
migration characteristics were affected by both the applied 
load and the rate of heat removal. The effects of the heat 
energy removal rate indicated that the concept of an energy 
balance should be considered in the development of a theory 


predicting moisture migration rates in freezing soils. 


apeniexh exArRye nage - ebay 796 boiz1BD mee 


asd ‘tain oF tuo BoLaaNe sou thie e 
sebve #412 noes fo eorsairsagazsts “soe texas 


{kA -ots2 erisasad foe feot failaga 40 anoessbaos 
al 


duo Peiaqss Siew efitoucieTxys atod pHok hns wIeF init 
eis+ 450de 2ure= lo ytrTOrer oc? 20d ted? havod as “ob B 

SeiSesut «ant poot iothe?q oF fees ot Bkwon gai 
Jou 695 od¥svasedO eigs Yo enottastati dt 
to nobse eh of 2p aenilabiep, aavseod booderehan — * 


pokisopia $30 SeiLor sua? atte? nist*do 68 


-cevip ots mime 


Sietetoe os teh+ pavers eenpery Ie thenZisqws ott 


: 
hecigqgn eS Aided yo belo lin Suaw coktaanataszedy agotse7pe is 

u 
+e@d Sd? Io 2Fo0I4e ac? .lovoeor cest Io oFey ORY Sen! Beek 


YPIOne oh Io Iggonon oft tod) Lotenipat stea Laveren ym 


= ; 


Ytoas? 6 tO tnsAqclovel of3 wl iistahbienoo ad bloods 
«elton palsenzi ni aeiai noiseapin sadetoe pat 39 ho 


a | 


a 


ACKNOWLEDGEMENTS 


The author would like to thank Dr. N.R.Morgenstern for 
suggesting the thesis topic and for the invaluable help and 


guidance given during all stages of this study. 


The author also thanks Messrs. 0O.Wood, A~Muir and 
J-Cyre for their assistance in both the experimental progran 


and the thesis preparation. 


Gratitude is also extended to to the many graduate 
collegues especially Messrs. R.W.Tweedie and J.S.Weaver for 


their assistance during apparatus set up. 


The author also wishes to thank Ms. S.Kowal and Ms. 
S.Searle for their help in the final preparation of this 


thesis. 


Finally the author wishes to thank Messrs. K.D.Flock 
and L.V.Medeiros for their help with the computing work 


carried out during the the experimental progran. 


vi 


ro? ‘ Q Amel pe 
is sie afdenbawas add 163 Ore siged © ate 
Naito 2id? Ve aonera tts vain os 
+44 aan 


tie ‘slowed .hooWs0 .ere20" ekasta “Cetb Pern 


rT 
: oof hl mn 4 
wbupord tetaeetsoqus os dfol 0) sonstalens aisdd 30% #12y- 


-as Ss ay oagy — ail 


’ + on 
efegbnay Yass et? of oF Safvetie cele al. wbys i thse | 


me 


40% Foveov.2... bus alboowT v.82 .Atee0R qifsiaages ete 
rar ie 


a toe andsosgae pilawh epnsielaee aes 
‘ev pe? 


obe has lawet.e <2 (api? o* eelletw Opis Godaun sat 


eid’ 16 NOdDRueqQoSY, ientt edt-as qlent Wsns to et 3% 
i ° ] 
’ ' a . J 3 - 


og 57.0.7 s#eneen Nass oF osteee, zon tre add yileaes . _ 
wzoe “palisqeoo ens Ajiw (gis sited 204 Setietheat bom 
| _ a 


—MbIPOIG (sIVSRiTegAe Skt en tocsvh suo ‘he 


TASLE OF CONTENTS 


PAGE 
Abstract esaenwnerteesvneeven ese ean es @Ovneeveneresne2e8 esse ee seeseneee#wwert#eneeets V 
Acknowledgements aeeetane eee eee seeowwvweesweeneesesesneveeeeveaene @Beeese vi 


Table of Contents seeevwvwrteznstwneveveeesvsseeeevseees eo 2 


aha Shas evita totaal ater Tal 


eeeeesevnee ses ix, 


Pr tOL eT AO LES ctarttag aw wie etatetadrats wie c-atate akeiatc chete 

DrStmOUer Ee GUILES Wyck. ac a see eke Oe Sie ee Sos sore cla, GX 

List of Photographic Pilla Wes aitite Ss ehars as a 7a ahete ia A tidratel Wei att Co: 

CHAPTER 1 PNT RODUGCT IGN sonst atst s)<\s/ eis ere ater erct als satel dee wieia sts 1 
si Permafrost and Associated 


Geotechical Problems (sf ac. as 00 en 6 1 


Ts 2 SCOpe sO SE Ret Phe Si Sele eats stele gists.» 
CHAPTER 2 LITERATURE REVIEW ac sccassec stessescsncss 4 
2a ted Freezing Behaviour cscccscersersee 4 
Carlo Mechanism of Frost Heaving ...... 6 
Ze te3 Water Pressures During Freezing . 13 


25165 Effects of Heat Energy Removal 


Rate =asaeveese@eenvpseeeeseee#eeesestsse ese &@ @ 19 
ig at Magnitude of sHeavetss shale ssc ces 120 
geen Basic concepts of Frost Heaving . 23 


HAS hy Ps Energy Budget saaesscneceasecencee 23 


A Rote TO Potential Gladvent.. oi «sss ee else oO 


CHAPTER 3 TEST PROCEDURE AND DESCRIPTION OF EQUIPMENT 
Bo oet MAtCLTIALS cece wus ssw eles sewevicces £0 
See Testing Program and Procedure ... 28 
34.5 Description of Equipment ........ 32 


vis 


uae ae 


: v. eee ee +4 Rhee tie 


Se Sew we ¥ cm 
fin ® te Ge a ce a ee 


ket yidthcermeeaaan eos 


» 


PIV oe ct wavs aia s iowa @eoe 


WE dedde duct usba eps cobiedes sate et ateerp hee MORRE Ben 

| Sewanee OOss Ke bi orth cee aw eenene men anrmeet 10:38 
chweeds veh BOMEIA deiqntpoFeds ei 

KoTtotdoart: | F aurea - 


tx We eee ew ee wreoeorrev €- 


OF £l we aiihs Selassie we nena <ks 

bereioaursA iar feo 7 LAs bat : 
t. . deawwesecencs Ensidcog? Leorinataste ; 
O° . dedbecwcenss ShRotTrpet7 Jo etene Cue 


B® (vite s ceaee eas Decbussevns WOTFSH SHU TARTS 
b eeew see a eee mee CIC) Vee prado. rt se | ie a 
WH Gavees bake was JE097 TO Meakabossn ~ Yo ag * (ees 
Ef « polsoory parca NATTA Cass 18 teW bal 95 
AP snavesees BAO) Poslagn % wioeIks DFS 
Esvomet yeuent tee oO = 901 TE cub st 
OD aw ee ek cod tps teense se oe see oes HOOT 
DS we wksawnse eyes S¥AGH VO Shut inpsr a0 
| ts . gaives! Fabet im 2iaonos SLenG <8 
é 

TB vnc tvewsdacdsee- nk sophie Yevont f.8ue 
| (BE awed tune deaees die hhexe (eitnegos Est.6 

re manent worraranes a, ate a@vag 0a Tea) t 
ere “es mies: ROE RIE Ree ee tye. 


t caipobe . 


CHAPTER PAGE 


CHAPTER 4 AE) fade Uda) Scie we, wate ge ere Si whelnie Seale ra stele watereietet) Ue 
4.4 BOESe. POOL aM feratec) stele ate wel siaia elses a oe 
452 RESULTS Ol SEf168 NDS sdvawcawmesana 49 
4.3 Results ~oh VSecies WNDS-B esse woeme os 27 
44 RESULTS SOL sSCG Iles aNDS=C esajewisieie oc OC 
4.5 Results of Sen LeS NDS =D seo oc ciate o OZ 
4.6 Conclusions From Observations ... 63 
4.7 Limitations of Experimental Data 71 

CHAPTER 5 ANA LYG2 3S OR TES TRDA TA Biterc ooo tere sare atatclcrese G04 
eA Temperature Data AnalySisS ....... 74 
Sore Pore Pressures Developed at the 


Frost Front eaeaas#2sesteeceesseesesaeaenaee ese & 84 


CHAPTER 6 CONCLUSIONS AND) RECOMMENDATIONS: <..24 0 cece GO 
6.1 Conclusions eseeevnww#eteeanaeneseoes#seeseaeaees 83 #8 8 @ 83 
6.2 Equipment Performance and 


PROCOGCULE sare sre wieielaieiwis ets oe ois) a slaceie C7 

6.3 Recommendations cscrcccsesssocses G2 

Vist OF ReCLCLeNCES San sc wee ow cas ome culm Seeiclssisn se sm awiew ce 94 
APPENDIX A PKOERIMENTAL- RESULTS 2e so atsiuls sss wslesiee sms 000)! 


APPENDIX B PROGRAM LISTING AND SAMPLE OUTPUT ........133 


viii 


7 


‘te anpbpeuchee anh a5i2a0 Ie. wature 


Co ripen Vaan eee hii 


We {iaees douse GRE ealgae 20 astuas x ee _S 
DO <5 os Hos D-H | 20l 192 1a 2s fpewi oe 
BA sacese eee G-cMh Zhtta? Jo -2tLuze? | Pit : - 
fe obe Bnoksevised® cost anolep inde? Dat a 

0 (een Lesaemiceuxa Jo. cooitsstaks rw 


area THet 30. 27Tevante 


UT 6 Salsa be ee ht ce bes 

UN cidecdus BFAVisns ote san sate rset i .2 : 
ei? Je begofsrst 26ers os3 S459 Se 

aba? “Fors . ¥ 


rt es eaed@eaetwaier ee ~* © oes 


BP ede ssenae GHOLTAQMS SEOOTR dhK FO DEM ORD 


Fe ols emt dee ten O40, 22 © ennrenaisvad fa! oa _ 
bis Pon ANIOTI4% jnana rays bed 


ee ehdm <w hae bare eet ae ads siphe sets 

CP wba ee eer ke wee enh debdaecdos f td 
OOo Swe ca wre ce hs bv uence Cae ned wbduthduldce BODRSSRROe to 3 
fithovceedecacsocebscsse, @tglash Lerwanete age A , 
Of Tscceeee THIMUO A2siNAY GMA OHTTZTA MARDONA 


List Of sfABLES 


TABLE PAGE 
Zs! SUG Olt, BressuL es cOi, 15 tbe ssmreaiem «5.0 8/4.9,5,5,0,06- 110 
Sia! PEOpebt1LES. 0.2 4D.6ViOls Sal te Rio tanrciem 0 ete tebeieieie se 20 
4. Resid GSlOF a5 Oise se aN Gru syelabe mses sinie gees, eieieieiatenerets ase. oF | 
ee. ReSULES, Of series (NDS=2, Ch. s cienieisiisiciein' sc siete = OO 
Ko ROESULES MOT oer VES GND 5-6 pes icleetelats ciel cllcleie ares sane on! 
4.4 Results Of Series NUS=D sios~ een wells sees se 0)! 
Ss tl Suction deyeloped at the Frost Front ......... 84 


ix 


DIST OY FIGueses 


FIGURE PAGE 


2.1 Theoretical Relationship Between Soil Water 
voLentiral and (rreezing IPoint depression w..c0<00.6 15 


2.2 Schematic Representation of Relationship between 
Heaving Rate, Heaving pressure and Particle Size 
RES ECTNGa tr OM 1Ce (bens Grow lima ths arate a ale eaclalees oriea eed 


2-3 Relationship between Heaving Pressures and 
Average Pptateter’ Glo GaCelons™ 4-05 ses sia seem eGl et 


Bale Gran, Size" Curve (Of Devon Sat wows wee smile Siaivieieaen 
3e2 schematic diagram of Experimental Set-Up ~..<ise<e0-. 33 
32) 00S.) Drawings. for Experimental Apparatus c.s666. <4 
Soy OMe eC a NGS aN VCa WO, tis clei sas ateie ales sials cists slale im open be 
21 Sauple Experimental Results NDS 7 .ecsescvsscccess 47 
ae OLO eka GO SVeTSUS A DD1 Veda GU oits.. sisiew ele als te slw syne 
4.3 Thermal Data Withee thes pone eee eS oe oO em aS Se 
4.4 Thermal Data ARS a BUSI PS 6 is SoA On 9 SO SO eee 
USO mR el est ReSultS NDS (6) 26's cies wim «sles ele sis cisa sels cle ciswiete 0 
BO PE Test Results NDS 60 <2. snc scence cece steele ance essa, 0D 
4.7 to 4.15 Total Heave and Segregational Heave Data .. 65 
4.16 Efficiency Verus Square Root of Step Temperature . 72 
5.1 to 5.4 Temperature Data AnalySiS s.ssessesssecesees 18 
5.5 Suction Developed at the Frost Front ..csseeeeeee- BU 


5.60 Schematic Diagram of Pore Pressure 
Histribution in the Soil Sample .....<.sceescsces- 86 


Oly ee ° 4 


as : ’ , yy : d 7 
aaugart A0veta 


gerne Lioe. udawtoy usdanoitAte? A bod: Snape 
Thad pk Zae eenne pokeensged $7209 yarseeat eee ot 


rs ciwee b> wee OP PO eee ees seeor em ta e352 207 


tne. § #etes sot privsed feapried qideao rte, 
enhorinen tt oO iw tombs eye 


PS *neeefeeesetrre ere He 


OF a@e eG ew Beeches ere +ire your oF to seid tng r? ahaa 


te? Eetaarraeqee 12, ah peed 9b tsagdoe 


ay eeeetnunee = i- 


df occ dus! HOST ECGA Teduoniisca? 209 epasbedd Th eae 


ai 
Oe ee bean ok ae 6 eo ‘soibat spnpdd oni toy f° 


o! 


ve ee eese @eaCenvnreve-*# , eu 2s lhe ey (5t*uSmi Ie sss aigape Va ad 


vi 


Oe ee ee ee ee ee ee Eao.] heztluy & eurazey oliys# bkOY ‘Se 


. 


- 
i ae Gee ee ee eee ec tie ee we ee 8 v ra Ted PAM 5204 tein Sect & 7 


Poa tie 


GP Sele cw dbecdedeuwcd beve Uf A0n-F OW V coee Peetene + 
ne Ie Fee eee oe ee latin aid aad ~wo. © BE aélozeh saat ioe 
BO Gaia Mee a elen rt ble de awae > Oe bw ae 2 ca S3tukes j2zeT 
Go &. SISO 4VFAH Tedoiteptape? bos sveeh fast er.e ay A 
Gv 4 S3¢e20qe6T ware to too8 Inu pe eyaeh youn o 2 NS ate 
BY ds pecascvesawes... ehdylenk stko orate 2eqEe® v.2 ott 
OB \odsssntasdass JSNOTT S209 sav Se heqoteved Maisie? < 


82025414, Sto) IW. CSIGEAT 6 2 a : 
an ee es slamee rue ‘eas nt netiealate oT, ’ i - 


LISTVOr PRQRUGRAPAIC PLATES 


PLATE PAGE 


1 Bott zanta |} Snadvysrtacat. ce eenses ta.s eee wc 56 


X1 


aoe@enrweeete eee 


7 ee 
2yradt: apnea 70K 30 
7 i 7 : a, ‘i Ain i 
ware Won 
enateensol 


- span) ed. tT evbigey bne ; 
- er DAS 2 


=| ; 7 


7 


2S 


- 
wm 


o. eae 


_ 


~o 


rs _* 


a | 


eelnmec 


41.1 Permafrost and Associated Geotechnical Problems 


With the rapid development of the industries associated 
with the recovery of natural resources found extensively in 
the Canadian North, an understanding of the behaviour of 
frozen, freezing and thawing ground is necessary to permit 
the safe design of the many and varied structures which this 


expansion demands. 


/Recent research has widened our understanding of the 
mechanics of cold region geotechnics enabling the design 
engineer to predict the behaviour of soils, frozen and 
thawing, accurately both obs) terms of strength and 
deformation. For problems involving thawed soil the 
conventional principles of temperate zone soil “Pence: are 
applicable. Work by Nixon and Morgenstern (1973) encompassed 
the problems of thawing ground and the development of a 
method of determining the magnitude of pore pressures 
generated during thaw and the prediction of thaw settlements 
was presented. This was a comprehensive theory which 
considered the two controlling dynamic factors of pore 


pressure dissipation and a moving thaw boundary. 


The characteristics of thawed soil were investigated by 


Roggensack (1976) and its behaviour is now well understood. 
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Freezing problems however are not as well comprehended. 
ohe petfects: of freezing a frost susceptible soil are nost 
clearly exhibited in a heaving of the ground _Surface. The 
magnitude of this heaving varies as a function of several 
factors such as applied load, stress history and grain size 
distribution. Further investigation of the soil after 
freezing reveals the development of discrete ice lenses 
which may grow to appreciable thicknesses. A decription of 
these lenses was given by MacKay (1972) in a study of ground 


ice occurrence in the high Arctic. 


The build-up of discrete ice layers 1s caused by the 
migration of water in the soil under the influence of a 
potential gradient created at or near the freezing 
interface. The prediction of the magnitude of this potential 
gradient has so far been impossible due to the complexity of 


the problen. 


If the heaving of the ground surface is inhibited, 
substantial heaving pressures may be developed, OT 
conversely if the overburden stress at the frost plane is 


high enough the growth of segregated ice is impeded. 


Another problem associated with the build up of ice 
layers is the effects of any subsequent thaw. The localized 
increase in water content of the soil implies a decrease of 
strength on thawing if drainage is hinderd by some factor 


such as low permeability of the soil. 
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To date no comprehensive understanding of frost heave 
behayiour is known and consequently the effects of the many 
controlling variables in the problem must be investigated. 
The purpose of this study was to examine the freezing 
behaviour of a frost susceptible soil under varying 
conditions of applied load and temperature distribution, and 


to correlate behaviour with those parameters. 
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CHAPTER 2 


LITERATURE REVIEW 


2.1.1 Freezing Behaviour 
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it is known .that most saturated soils under low 
stresses exhibit an affinity for water during freezing and 
the formation of lenses and layers of ice occurs. Coarse 
grained soils may expel water on freezing if drainage is 
permitted, (MacKay 1971, McRoberts 1972). In reviewing the 
geotechnical properties of freezing ground, Anderson and 
Morgenstern (1973) stated that ice lensing results when a 
SsOil has the ability to supply water to an active ice front 
for a sufficiently long time to grow an ice lens. Frost 
action is contingent on the existence of certain factors. 
These are a frost susceptible soil, an adequate moisture 
supply and a sufficiently low temperature to cause some of 


the water to freeze. 


The effects of frost action are displayed mainly in the 
form of vertical displacement of the soil surface. eee 
problems also occur, @.g9.,; the increased soil moisture 
content due to water migration during the freezing process 
is reflected in a decrease in soil strength after thaw 
occurs. Vertical displacement of subsurface objects such as 


boulders, piles and posts are also of major concern. 
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volumetric increase of a soil sample on freezing is due to 
two components: 
1) The volumetric increase of in situ pore water due to its 


phase change. This may be expressed mathematically as 


dh dx 
—— = 0.09n —- Zon 
dt dt 

where Nn = porosity of the soil 


x = depth of frost penetration 

t = time 

h = heave due to freezing of insitu per water. 
2) The volumetric inqrease due to moisture transfer and its 


subsequent phase change. This is given by 


dh 
—- = 1.09V Zoe 
dt 
Where h = heave due to the formation of segregated 
ice 
v = velocity of the water arriving at the 


frost front due to a potential gradient. 


It is the latter which is of major importance in 
engineering considerations since the suppression of heaving 
due to freezing of insitu pore water implies the application 
of loads high enough to inhibit crystallization completely 
and these magnitudes are well beyond those of engineering 


interest. 


Taber (1930) first showed the phenomenon of moisture 


migration in an experiment in which a soil with a pore fluid 
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of benzene was frozen. The sample heaved on freezing and yet 
benzene itself exhibits no volumetric expansion on freezing. 
He concluded that the build-up of frozen benzene lenses and 
the resulting volume increase and heaving pressures were due 


+o moisture transfer. 


The affinity of a freezing soil for water is due to 
soil water suctions set up at or near the ice/water 
interface and it is these negative pore pressures that are 
of great interest. Theories have been postulated within the 
framework of capillary physics and thermodynamics, both of 
which give an insight into the mechanisms which cause these 
Suctions and the resulting water movements. However none as 
yet completely describe the ice/water interface phenomena 


adequately. 


2.1.2 Mechanism of Frost Heaving 


Whenever a negative temperature gradient is imposed on 
a soil the volumetric heat is removed causing a drop in bulk 
temperature. When all the volumetric heat is removed and the 
bulk temperature is 0°C, nucleation of ice occurs with a 
subsequent release of latent heat. If supercooling of the 
pore water occurs nucleation will occur at some temperature 
lower than O0°C. The magnitude of the supercooling effect 
depends on several factors. Available evidence indicates 
that 2°C to 8°C of supercooling may be required in the 


laboratory but little or none in natural Surroundings. In 
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coarse-textured soils the freezing interface closely follows 
the 0°C isotherm in its downward course. In fine textured 
soils it may be as much as several degrees behind due to 


freezing point depression (Anderson and Morgenstern 1973) 


It is now well established that soil particles are 
surrounded by an absorbed layer of water with different 
properties from that of the pore water (Anderson and 
Morgenstern, 1973). The thickness of this film is a function 
of temperature (Anderson, 1968). This water which does not 
freeze in a similar manner to the pore water is known as the 
unfrozen water content. To date the unfrozen water contents 
have been correlated to surface area, (Anderson, Tice and 
McKim, 1973), Atterberg limits, freezing point depression, 
Clay mineralogy and activity ratio, (Dillon and Andersland 


1966.) 


The equilibrium of the unfrozen water content may be 
described by the Clausius-Clapeyron equation applicable to 


this situation 


Dp =bH es 
ws; t 2e0 
Dt TAV 
Where Dp = change in pressure across the freezing 
EEOnt 
Dt = change in temperature across the freezing 
front 


AH, = The diminished latent heat of freezing 


soil water 
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The freezing temperature of the soil 
Ov = change in specific volume during phase 


Change transformation. 


The factors determing the solidy/ liquid boundary on the 
phase diagram for ice in frozen soil involve the magnitude 
of the reduction in the freezing temperature T and the 
diminished latent heat of freezing soil water, Ales - Data 
are not available to fully resolve the question, but 
estimates indicate that of the two opposing effects the 
decrease in Ales predominates, (Anderson and Morgerstern 


1973) 


The application of pressure tends to increase the 
thickness of the unfrozen water layer. Evidence that the 
thickness of the unfrozen water layer does increase with 
increasing pressure is provided in the observation of 
Hoekstra and Keune (1967) that the electrical conductance of 
a frozen clay water mixture increases dramatically with 


pressure. 


As more heat is removed and after ice nucleation the 
freezing interface descends. As this occurs one of two 
things may happen. In one case the pore water may freeze in 
situ aS growing ice crystals invade, fill the soii pores and 
engulf the soil grains. In this case the development of 
heaving forces results entirely from the expansion of in- 
situ water due to phase change. This is the mode of 


behaviour in the situation of freezing a coarse grained soil 


| oft fo yashatod Kivu pit Gee ait) eer 


Shutiapea adtesviornt [ion pesert'* ab epi ae cones 
eth fos 7 aap te794 09 phiseoit eds ok me 
By60 . ai qaataw eine poixeg32 to Ssed Saaee 
tod )<neiyaetp “sit ev loees (ilut .o* a1dukeave: 
eds aysetis pminesgh ows oft Go edt stsoibad 28 
wiev2tepiah fest Aces dnA, 2° ten tagboard (ye ea 


ait aesetont o7 .aheoo? salaseIg, 26 cased tone fst 
of? seus gnuppu v5 raat kt nenenantrss Hy area 
fiiu gasostonr 2 ‘ba Ear so ta penon tay edd. ae 
Yo nol#etasado BAe. nv) pebivorg ei sm23e3g Bee : 
30 sca6suubags fsointoslio edt, 509 icaen} enue? Bas 


Astin NI Léodtsrsrh Sie no aie a7 sxi a) Peet ytis 


edi ooktestonn “ang 793%6 haw heveae3 Ri Teed s708 * 
out Io S30 sTunD0> etd} eA ae, yon Piesis 


f 
ig 


mae 69% gaoq' Sis sanoon 


% 


in a closed system configuration. Alternatively, segregated 
ice lenses may form in a plane of favorable temperature near 
but a little behind the zero degree isotherm. The ice lenses 
may thus enlarge by the addition of ice formed from upward 
migrating soil water. It is this upward flow of water and 
the potential and thermal gradients which control this 
Migration that are of special interest. The frost front 
moves in accordance with the relative balance or unbalance 
between heat energy brought to and removed from the frost 
front. An advancing frost front results in the engulfment of 
soil particles by ice crystals or in situ freezing. If the 
rate of heat removal equals the rate of heat supplied, the 


frost front becomes immobile, i.e. a stationary frost front. 


Once the freezing front becomes stationary and water in 
sufficient quantities is available to the front,the ice 
crystals grow by freezing the adjacent water. Due to _ the 
rigidity and perfection of the ice crystals, soil particles 
are rejected by the ice, resulting in the segregation of ice 


and soil particles (Anderson, 1968). 


The heat balance at the frost front is a function of 
several modes of heat transport. Conduction of heat through 
the frozen zone is the sink of heat removal. Conduction of 
heat towards the frost front through the unfrozen zone and 
heat released by the porewater and by any water flowing into 
the sample are the major sources of heat energy. Heat 


released by the porewater is volumetric heat and latent heat 
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of fusion. The latent heat of fusion is the major source of 


heat (Nixon 1973). 


There are two commonly adopted explanations for the 
mechanism of ice lens growth. Consideration of capillary 
physics, where the soil pore is analogous to a fine 


capillary in which ice is nucleating yields the Kelvin 


eq uation. 


P; mie oe 260i, 2.4 
iw 
Where P, = pressure on the ice, overburden pressure 
Pw = pressure in the porewater 


O., = pressure differential across the ice/water 
interface 
C.=-tadius of the ice/water interface. 
Many forms of this equation have been developed for various 
configurations e.g. unsaturated soils where the effects of 
air tensions must be considered. The formulation above is 
the basic equation for the pressure difference across an 


ice/ water interface for small ice crystals in their own 


melt. 


Penner (1959) postulated a complete capillary model 
relating pore size and associated suctions. Assuming the 
size of the soil pore determines the size of the advancing 
ice crystal, the quantity fc, in equation 2.4 may be 


interpreted as the radius of the pore. A part of the soil 


pore is occupied by absorbed water which probably places an 
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additional size limitation on the ice crystal. As the frost 
front progresses, its propagation may be temporarily halted 
if the radius of the ice crystal is greater than that of the 
soil pore. However, heat is still being removed from the 
system and since the ice crystal cannot enter the pore space 
supercooling of the pore fluid ahead of the frost front will 


occur. The depression of the freezing point is given by the 


relation 
Sama eae ermal ear aeea ar Phe ey 
Cae 
Where T = freezing point in degrees Kelvin 


T,= normal freezing point 


Ve= specific volume of water at the freezing 


temperature 

6.= stress difference across the ice/water 
interface 

[y= radius of the ice water interface 

L = latent heat of fusion of water. 


AS the porewater is further cooled it may be noted from 
equation 2.5 that the radius of the ice water interface must 
decrease to maintain equality of both sides. Supercooling 
will continue until the radius has become small enough to 
enter the soil pore. As the radius of curvature of the ice 
crystal diminishes, the Kelvin equation (equation 2.4) 
predicts that PB-P, will increase, and since P, is 
effectively the overburden pressure and hence not a variant, 


ed must decrease. This decrease in the pore water pressure 
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will generate the suction gradient which draws water to the 
freezing front. The freezing of this water results in the 
growth of ice lenses. This theory has several shortcomings 
and is not capable of describing all freezing conditions. 
Some of these are: 1) The capillary theory is based on a 
uniform pore size or particle size distribution. In nature, 
virtually all soils exhibit a non-uniform particle size 
distribution and consequently difficulties arise in 
asSigning a typical radius for use in the theory, 2) The 
capillary theory is a static theory based on the premise 
that pore radius or grain size distributicn is the ee 
parameter controlling the ice lensing capability of the 
soil. The theory does not account for any other factors 
affecting the freezing soil system e.g. temperature gradient 
or stress history and is incapable of accounting for any 


changes in these factors. 


An alternative mechanism for moisture migration is 
Known as the thermodynamic model. As discussed previously a 
soil particle has an equilibrum value of thickness of 
adsorbed water layer for a given set of conditions of 
temperature and pressure. Considering a soil particle at the 
frost front, some of the adsorbed water on the cold side of 
the particle may crystallize and leave the film. In order to 
maintain the equilibrium film thickness water molecules fron 
the unfrozen void water are absorbed by the film. A process 
is thereby set up to move water molecules to the freezing 


front. This process would also imply the exclusion of the 
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soil particles from the ice structure giving lenses of pure 


ice on freezing, as found in nature. 


Knowledge of the pore water pressure distribution 
through the soil sample is of utmost importance in the 
solution of the frost heave problem. In the preceding 
chapter two mechanisims were discussed whereby a potential 
gradient may be set up on freezing. Investigation of this 
potential gradient by several authors has thrown some light 
on the problem. Williams (1967) in a series of experiments 
in which the pore pressure was measured as the frost front 
penetrated through the sample,showed that a negative pore 
pressure was developed at the frost front and concluded that 
this was responsible for the potential gradient causing 


moisture migration into the sample. 


Calculations using the Kelvin equation yield values of 
pore press developed at the frost front and these are found 
to vary inversly with pore size (Penner 1958). Williams 
(1967) presented data showing that the value of pore radius 
obtained from air intrusion experiments could be used to 


predict pore pressures developed at the freezing front. 


Sutherland and Gaskin (1973) showed that pore pressure 
at the freezing front could be reasonably predicted using 
grain size and capillary pore size, and their work lent 


support to Everett's (1965) proposal that in calculations 
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the fy value to be used is associated with the effective 
size D, rather than the mean particle size Dg. However, 
the authors concluded that predicted values of the drcp in 
porewater pressures using values of pore size obtained by 
the air entry method grossly underestimated the actual 


measured values. 


In the frozen zone the movement of water which occurs 
under the influence of suctions in the frozen zone is 
thought to be induced by the different energy levels of the 
unfrozen water and ice in the soil matrix. Harlan (1974) 
equa‘ed the Gibbs free energy difference between the two 
phases with the potential of the water in the soil in 
equilibrium with ice at a constant vapour pressure, 
temperature and pressure. Provided the ice crystals are not 
too small and have properties not unlike those of bulk ice 
at atmospheric pressure, the vapour pressure of the ice is 
independent of the amount of ice in the system and similarly 
the unfrozen water content is independent of total water 
contents. The development gives a relationship between 
temperature below 0°C i.e. freezing point depression and 


soil water potential as in Fig 2.1. 


Harlan predicts a suction profile in the unfrozen water 
in the frozen zone and perhaps continuing into the unfrozen 
zone. The flow of moisture in the frozen zone has been 
demonstrated and documented by Hoekstra (1966). Due to the 


very low permeability of frozen soil the amount of flow 


jaeeeor er, ae he ' ey? 76g 
wt gob eid a0 eater. ee ; 


qa fenisido ssie a6 io ‘aie 
Teuton $45 hate ldeedebrit yleaore beau had mee tf 

4 ; 
‘tame vog) Hao» saan aexora att Ath ‘2 
jqpae 40 esnoodser ost web am 


ont Qo afovel yprens ts 297R Lh oi? vd ‘py dvbot ad o* 


ce 


agn0n0 “otie satiew 3 


#2. Sros nayvyozl. sé ni 


(prety) nslvsH ,-xiaioe ‘Line edd we aad bos soteu # 
ov+ sdf Hasuded /gornsis7 20 yorsu8 pent 2dGh0 eat bes : 
at . fzem) sft wi 19708 8d] lo | Cativetey oft atneir mt 


eatezeiq 2u0qRV Fi eIaaeS 5 45 ok doi rab 
fon 975 eid s2429, ied ot, bor lypaa’’ tie siveniph > bee # 
ear Aid a0) shorts erin son, B6fIt9CoTs eyed fee ae 
ak) sok. edt to ‘emeent9 rb OF GY ‘pay ewesone seen i, 
yireiisa ia bas Stays 44 nt BDk 20 +nuswh slg to ini | ak 
sovev [6907 920 fasbasyehdi mk Sietoos aste¥ eric a 
iAewi od qidenotts| 43 6 asvip | ane: lid SAY areas a8 
(bne noignexosb maieG putea 7? out 3° 0" vurett m ” pe 
1S vit ohne Leignstog 39. 


or perro treet 
Dinas par lana 
a ee “ss aie he HOSTS 


<a 


t) “nedeveoh 


at 
ee 


5 


-15.0 
oO 
° 
= 
= 
o 
%-10.0 
x 
a 
J 
(=) 
e 
= 
(e) 
a 
-5.0 
oO 
= 
N 
uJ 
uJ 
og 
Lo 
0.0 


-1x 10> -1x104 -1x10 -1x10® 
SOIL- WATER POTENTIAL - cm. of H,0 


-1x 102 


Fig. 2.1. Theoretical Relationship Between Soil 
Water Potential and Freezing Point Depression. (Harlan, 1974) 
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Table 2.1 Shut-off Pressures of Soils 
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induced by these tensions is negligible in engineering 


terms. However the mode of behaviour is of interest. 


By examination of the above, the concept of a discrete 
interface advancing through a soil comes into doubt since 
not only is there formation of ice at the 0°C isotherm but 
also at some distance behind it at a position of favourable 
heat conditions. This greatly complicates the problen, 
however it is believed that for engineering applications the 
post 00°C isotherm freezing may be neglected and the 
Simplifying assumption that all of the ice formation occurs 
at the interface be made. It is also assumed that any 


suctions generated originated at or near this interface. 


2.1.4 Effects of Applied Load 


ess 


The fact that increased load can depress the heave rate 
has been appreciated by workers in the field of frost action 
for some 40 years and there are several well documented case 


records which illustrate the phenomenon. 


When an ice/water interface advances through a freezing 
soil it can be observed that water is either attracted to or 
expelled from the freezing front. It is found that for the 
same soil type attraction occurs when the applied load is 
low and that expulsion conditions exist when the applied 
load is higher. Moreover, there exists a unique stress, for 
a given soil type and stress history at which no flow into 


or out of the sample will occur. This stress is called the 
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shut-off pressure. 


Beskow (1935) was the first to document that increased 
pressure causes a decrease in heaving rate ina silty soil. 
Beskow manipulated the stress level by either changing the 
surface load, the pore water pressure or combinations of 
both and found that the heaving rate was identical if the 
difference between applied stress and the _ pore water 
pressure remained constant. The larger the difference 


between the two stresses the slower was the heaving rate. 


Linell and Kaplar (1959) showed the existence of a 
Shut- off pressure for a variety of soils. Their results are 


presented in Table 2.1 


Ot her supporting laboratory evidence presented by 
Anderson (1973) is aiso summarized in Table 2.1. These tests 
showed that when the applied load was lower than the shut 
off pressure water was attracted to the freezing point and 
when the applied load exceeded the shut-off pressure, water 


was expelled. 


Theoretical justification of the existence of the shut 
off pressure may be gained from a consideration of the 
capillary model. Equation 2.4 predicts that a pressure 
difference exists across the curved surface of the ice water 


interface due to the energy of the interface. 


Two conditions are predicted by this equation. If the 


applied load on freezing is less than ( 26, ) #£=water 
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migration occurs towards the interface and an ice lens 
grows. If however,the stress difference is greater than 
( 20w ) then ice will propagate through the soil pores and 
water may be expelled. This stress ( 26, ) is, then, the 


shut off pressure. 


Arvidson (1973) equated (20 wy) with the shut-off 
pressure and attempted to predict the value of shut-off 
pressure using typical values of ox, and ¢,,,however the 
agreement between predicted and observed was not good and 
Showed the limitations of assuming values for parameters 


which cannot be accurately assessed. 


Takashi (1974) in a series of experiments in which the 
rate of penetration of the frost front was kept 
constant,showed the different modes of behaviour of moisture 
migration due to changes in applied load. Values of shut-off 
pressure obtained ranged from 1.2 Kg/cm? for a silt to 3.5 


Kg/cm= for a silty clay. 


An example of a case record in which both total heave 
rate is known and where the freezing rates are given is 
presented by Aitken (1974) for a field test program 
conducted near Fairbanks, Alaska. These tests were conducted 
by monitoring the total heave experienced by a series of 24 
foot square loaded platforms. Allowance for heave due _ to 
freezing of situ water in the total heave results gives the 
segregational heave rate with respect to applied load. 


Aitken showed that the segregational heave rate decreased 
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almost linearly with increasing applied load. 


The effects of the rate of removal of heat from the 
frost front are uncertain. Beskow (1935) demonstrated that 
the heaving rate was not always influenced by the rate of 
frost penetration. In his experiments the temperature was 
varied between -3°C and -10°C above the sample without 
affecting the heave rate. The colder temperatures, although 
increasing the rate of frost penetration, did not increase 
the heaving rate. Penner (1959) however, showed 
experimentally that the rate of heave was a function of the 
rate of heat removal using cold side temperatures of -2°C to 
-6°C. YThe difficulty in assessing the overall effects of 
heat and moisture flow arise not only from the complexity of 
the mathematics but also from the lack of experimental 
measurements of heave rates, heat flow, moisture flow, 
temperature distributions and moisture tensions while ice 
lensing is in progress. In the absence of such information, 
the quantitative treatment of the combined heat and moisture 


flow appears very difficult. 


Tagaki (1970) presented an analysis of ice lens 
formation. This work formulates in great detail the heat and 
moisture flow equations in a layer of unfrozen soil overlain 
by a layer of ice. This theory does not predict the 


cessation of the growth of this ice layer or the advance of 
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an in-situ frost front beyond its initial position. Although 
extremely rigorous in presentation, the theory in its 
present form is not capable of predicting the magnitude of 


frost heaving likely to be experienced in the field. 


2-1.6 Magnitude of Heave 


The magnitude of the heaving observed during freezing 
is due to the combination of heave due to freezing of in- 
Situ water and heave due to moisture transfer. Fora 
Stationary frost front the former of these is reduced to 
zero and hence only segregational heave occurs. With a 
stationary frost front the component affecting the magnitude 
of heave is a function of the magnitude of suction generated 
at the frost front. If the suction is decreased by some 
means the rate of water flow into the sample would be also 
be decreased hence lowering the magnitude of the rate of 


heave. 


The effects of applied load and heat flux removal for a 
given soil have been reviewed previously, however the 
effects of a different soil type for conditions of constant 


applied load and heat flux removal must be considered. 


The most obvious factor is that of particle size. 
Penner (1967b) summarized heaving behaviour in a _ schematic 
graph (Fig 2.2). In zone 1 the heaving behaviour is subdued 
since at this grain size range low values of suction would 


be generated and although the permeability of the material 
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MAXIMUM HEAVING PRESSUR 
CURVE OR SUCTION INDUCED 
CONTROLLED BY PORE 
AND PARTICLE SIZE “ 


HEAVE RATE OR 
MAXIMUM HEAVING PRESSURE 


Coarse grained Medium Fine grained 
TEXTURE OF SOILS 


Fig. 2.2 Schematic Representation of Relationship 
Between Heaving Rate, Heaving Pressure and Particle 
Size Resulting from Ice Lens Growth. (Penner, 1967) 
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Fig. 2.3 Relationship Between Heaving Pressures 
and Average Diameter of Fractions. (Penner, 1967) 
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would permit high flow rates, the potential gradient would 
not he sufficient to draw in an appreciable amount of water. 
zone 3, however, is the opposite case; i.e. high suctions 
would be generated but due to the combined effects of 
cavitation and low permeability water flow would be 
restricted. Zone 2 is the optimum for maximum heave due to 
ice segregation. A combination of moderate suctions and 
permeability render this particle size very frost 


susceptible. 


The effects of grain size alone on the suction 
developed at the frost front was demonstrated very clearly 
by Penner (1967a) by the freezing of glass beads of varying 
diameters and measuring the associated suction gradient. A 
curve of load diameter versus water potential was obtained 


and is given in Fig 2.3. 


Another factor governing the rate of heave is the 
availability of water. As stated previously, ice lensing 
results when a soil has the ability to supply water to an 
active ice front for a sufficiently long time to grow an ice 
lens, (Anderson and Morgenstern 1973). If that ability to 
supply water is removed the freezing process essentially 
becomes a closed system and the heaving behaviour is greatly 
subdued. McGaw (1972) correlated heaving rates with depth to 
the water table i.e. moisture availability, on four 
different soil types. The results showed that if water was 


made virtually unavailable by lowering the water table 
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suffiently, heave could be virtually eliminated. 


2-2-1 Basic Concepts of Frost Heaving 


The development of a comprehensive theory of frost 
heaving nust embrace the effects of the external factors of 
applied load and heat flux removal on the moisture migration 
behaviour. Consequently we must investigate the heat energy 
balance at the frost front and attempt to relate it in some 


way to the rate of moisture migration into the sample. 


2 Energy Budget 


The concept of an energy budget is a valuable framework 
on which to base any development of a frost heave theory. 
Although there is conflicting evidence on the subject, it is 
felt intuitively that the rate of heat removal must play 
some part in the rate of moisture migration. A combination 
of Beskow's (1935) results and those of Penner (1972) ed 
Kaplar (1970), with the suction gradient increasing with 
heat flux removal up to a certain level, and then becoming 
constant at this magnitude may best describe the effects of 


rate of heat removal. 


The transport of energy takes place under the influence 
of five factors; conduction, convection, radiation , mass 
transport and latent heat of fusion of the migrated water. 
Convection and radiation may be immediately neglected as 


being of no importance to this problem (Harlan, 1973). The 
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effect of mass transport may also be omitted as the 
magnitude of heat energy involved in this term is very small 
Since the temperature of the water flowing into the sample 


is very close to 0°C. 


The heat energy balance for a stationary frost front 
may be expressed mathematically by: 
Heat out by conduction - heat in by conduction - 
latent heat of fusion 7 0 


The heat out by conduction is given by 


dt 
oe 2-6 
az, 
Where Koo = The thersal conductivity of frozen “soil 
T = temperature 
z,= distance in the axial direction the frozen 
zone 


Similarly the heat conducted to the frost front by the 


unfrozen soil is given by 


at 
Py = a Lit 
dZa 
Where Ky = the thermal conductivity of the unfrozen 
SOx1< 


Z,= distance in the axial direction in the 


unfrozen zone 


Assuming a Stationary frost front, the heat brought to 


the frost front as latent heat of fusion may be expressed as 
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Where L = latent heat of fusion 
Ve=FVEVOCITy Of "the Water arriving ~at the 


Trost rront. 


The velocity of the water assuming at the frost front 


may be given as 


v = ki ea 
Where k = the coefficient of permeability 
i = the hydraulic gradient through the soil. 


Hence the sum of the heat transfer terms at the’ frost 
front may be expressed as 
Q,-2.-A7 9 AG 
The ineguality Sign being included to account for any 
additional energy considerations which may arise due to 
other components. Enumerating the terms gives 
dt dt 
K, —-K,—-LKi 7 0 Zee 
obras za 
The hydraulic gradient, i, may be more rigorously 
defined as the gradient of the total potential, i.e. ® , 
thus 2.11 may be re-written as 
dtm aval ag 
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The rate of moisture migration is governed by the laws 
of soil water flow. Water flows due to a gradient in total 


potential. Total potential is defined as 


® = yze (yw ve -& ) Baa2 

Where 6 = total potential 
¥24 = position head or gravitational head 
wy = pressure head or pressure potential 
é€ = osmotic pressure potential 


= adsorption potential. 


In this problem the osmotic potential and adsorption 
potential components may be neglected reducing expression 
Date £0 

P= yee} 25,13 
However, the presence of tensions in the pore water tends to 
complicate the problem, hence some formalism in definition 
and terminology is necessary to avoid confusion. The status 
of water in a soil can be expressed in terms of free energy 
relative to free, pure water (Aitchison, Russam and 
Richards, 1965). In cases where the free energy is less than 
that of pure water under atmospheric pressure the terms 
suction and negative pore water pressure are often used. The 
total potential of soil water is the potential in pure water 
that will cause the same free energy at the same temperature 


aS in the soil water. 


An alternative definition of total potential is the 
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work done per unit quantity to transport reversibly and 
isothermally an infinitesimal amount of pure water from a 
pool at a specified elevation at atmospheric eebe free the 
point in the soil water under consideration. At this point 
the balance of work done by the pore water must be 
investigated. In the definition given above, the potential 
is expressed as the work done transporting water into the 
sample from a pool at a specified elevation, however, in the 
frost heave problem there are two components of work: 1)that 
work done due to the flow of water under the influence of a 
potential gradient and 2) the work done in lifting the 


frozen soil and the applied load. 


It would appear that the greater the work required to 
lift the applied load the less will be the energy available 
to transport water at a given temperature gradient. This in 
general terms appears to be the action of externally applied 
loads. Consequently an experimental study of water flux 
during freezing with varying applied loads and rates of 
freezing was undertaken to investigate ' the inter- 


relationship between these factors. 
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Testing was carried out uSing Devon Silt, the 
properties of which are summarized in Table 3.1 and EEG 35:35 
This material was chosen for its high frost susceptibility, 
availability and its similarity to silty materials found in 


northern regions. 


The samples were prepared as a slurry at a moisture 
content of about 50%{1 1/2 times its liquid limit) for ease 
of handling. To ensure saturation the slurry was subjected 
to vigorous vibration, and a vacuum of approximately 680 aon 
of mercury was applied to the sample. This proved to be a 
very succesful method of de-airing the slurry as B values of 
unity were obtained both before and after several 


consolidation increments had been applied. 


3.2 Testing Program and Procedure 


The objective of the program was to investigate the 
moisture Migration characteristics of Devon Silt with 
respect to stress level and rate of advance of the frost 
front. Loads were applied using a weights and hanger system 
and the different heat flux conditions were obtained by 


applying a constant step temperature to the sample. 


Once the slurry had been placed in the apparatus, and 
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DEVON SILT 


LIQUID LIMIT OSV 


PLASTIC LIMIT ——= 165) % 


SPECIFIC GRAVITY —— 2.66 


% PASSING 200 SIEVE— 87% 


De — 0.0002m.m 


Table 3.1 Properties of Devon Silt 
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assembly completed, the test equipment was set ip inva cold 
room with a mean air temperature of 1.79C. The sample was 
then consolidated to the desired applied load. Temperature 
equilibration was ensured by leaving the sample for a period 
of time (usually 18 hours) in a styrofoam cabinet inside the 
cold room while each consolidation increment was applied. 
The B test to check the degree of satuation, if desired, was 


carried out at this stage. 


Prior to freezing the soil sample, the desired freeze 
test boundary conditions were imposed. Bottom drainage was 
closed and an unrestrained displacement boundary condition 
applied. The sample had free access to an external water 
source through the top plate of the apparatus. The volume of 
water leaving or entering the soil was measured by a volume 


change indicator system, described below. 


To eliminate problems of friction and adfreeze the 
samples were frozen from the bottom upwards by imposing a 
negative step temperature on the base of the sample. This 
temperature was held constant throughout the test and the 
resultant freeze rate was a function of sample moisture 


content and magnitude of the step temperature. 


The temperature distribution and movement of the 09°C 
isotherm was monitered continuously. The magnitude of heave 
was also measured. 


-3 Description of Equipment 
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The experimental apparatus was a Permode (Permafrost 
Oedometer) as described by Roggensack (1976) with slight 
modification to permit the measurement of temperature in the 


soil sample. 


The permode was intended primarily for work with 
thawing soils, however it lends itself well to work on the 


freezing of soils. 


A schematic diagram of the equipment set-up is shown in 


Pig es es 


The permode (Figs 3.3 to Fig 3.7) is a 100 mm diameter, 
teflon lined cylinder. The outer jacket of the cell was 
machined from thick walled P.V.C. pipe thereby preventing 
lateral strain during the application of load. The load was 
applied by weights on a counterbalanced hanger, and was 
transfered to the soil by the top piston. The seal on this 
piston was a double o-ring seal with grooves for the rings 
machined in the piston to increase the efficiency of the 
seal and to reduce side friction. The seal proved effective 
up to pressures of 1.4 Kg/cm2 though during the testing 
program the loading increments were held well below this 


level. 


A coating of high vacuum grease on the teflon reduced 
friction between the membrance and the walls of the cylinder 
to almost zero. This was deduced from the results of the B 


test carried out at various stages of loading. The results 
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of the B tests carried out on Sample NDS-6 are included in 


Chapter 4. 


\ 


The bottom seal was an O-ring and groove construction, 
effective to pressures above 1.8 Kg/cm@ due to the thickness 
of the O-ring. However, as mentioned before the efficiency 
of ths seal was never tested to its limit during the program 


because of the low load increments that were applied. 


Drainage was provided at both top and bottom through 
sintered steel porous’ stones. This is an improvement over 
previous apparatus as time for complete pore pressure 


dissipation after load application is greatly reduced. 


The base plate acted as a heat sink to induce freezing. 
It contained the coolant circulation maze which was designed 
to provide as high a flow rate as possible without denying 


adequate distribution of the fiuid throughout the base. 


Thermocouples were placed in the side walls of the 
cylinder to measure the temperature distribution in the soil 


sample. This system will be described in detail below. 


Freezing from the bottom was adopted for several 
reasons. The problems of adfreeze and side friction and the 
resulting restraining forces induced by trying to heave a 
plug of frozen soil upwards, (as in freezing from the _ top) 
was considered to be a major problem. Jf freezing was 
induced from the bottom upwards the soil being heaved would 


be unfrozen and hence side friction would be minimal and 
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adfreeze non-existent. 


Freezing from the base also eliminates the problem of 
having to know, very accurately, the sample height since the 
distance from the plane of application of the step 
temperature to each thermocouple is the same in each 
experiment, and hence temperature/time data may be easily 
correlated to observe and compare rates of penetration of 


the freezing front. 


The negative step temperature imposed on the base of 
the sample was attained by the circulation of a coolant 
through a maze in the base plate . The fluid temperature was 
Maintained by a constant temperature bath to 0.1°C. The 
fluid used was an anti-freeze/water mix (50% by volume) 


which was pumped through the maze at a rate of 80ml/sec. 


Thermocouples set into the walls of the cylinder were 
used to monitor temperature. A Honeywell 24 channel strip 
chart recorder recorded temperature distribution through the 
sample with time. Accuracy of 0.3°C was possible with this 
arrangement. However a slight time lag between the soil 
temperature and that measured will occur due to the low 
thermal conductivity of the membrance and the grease. This 
time lag will be constant for each set of tests and is thus 


a negligible source of error. 


Six thermocouples were installed at 10mm intervals and 


were numbered 7 to 12 from the base wpwards. Four other 
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thermocouples were employed to monitor coolant bath 
temperature, coolant return temperature, room temperature 
\ 


and the temperature inside the styofoam cabinet. 


Heat flux removed from the sample may be calculated 
from the above data. The ae difference between the 
exit and entry lines of the heat exchanger is directly 
proportional to the heat flux or alternatively knowing the 


temperature gradient in the frozen zone, and the thermal 


conductivity of frozen soil, the heat flux may be calculated 


from 
dt 
O=Ke-—A Sel 
dz 
Ke- thermal conductivity of frozen soil 
A - sample cross-sectional area 
dt - temperature gradient 
az 
Pore water pressure measurements were recorded using a 
pressure transducer calibrated to 0.017 Kg/cm2. These 


measurements were taken during B tests to ensure 100% 


saturation and complete load transfer to the sample. 


The volume change of the sample in terms of pore fluid 
entering or leaving the sample was measured in a 10cc volume 


change indicator as in Fig 3.8. 


Vertical displacements both during consolidation and 


freezing were monitored using a 6 volt excitation LVDT 
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(Linear Voltage Displacement Transdocer) accurate to 
0.025mm. The output was recorded on a Hewlett Packard data 


aquisition system employing a digital voltmeter. 


The apparatus, although not originally intended for 
this application, performed extremely well. Side friction 
was found to be negligible and no other major problems were 
encountered during testing. However, slight difficulty was 
found during set up. Since the slurry when first placed in 
the cylinder has zero strength the sample was unable to 
stand unsupported when the permode cylinder jackets were 
renoved to install the top piston and O-ring seals. This was 
overcome with the aid of an assistant who would support the 
slurry and membrane while the outer jacket was removed. The 
O-rings were then placed and sealed and the cylinder was 
replaced and tightened to complete the operation. Distilled, 
de-aired water was then circulated through the drainage 
system to remove any air bubbles which may have gathered 


during set-up. 


Overall performance was good. The only delays were 
caused by the failure of ancillary equipment such as_ the 
temperature bath, data aquisition system or failure of the 


room refrigeration system. 
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4,1 Test Program 


The aim of the test program was to obtain the moisture 
migration characteristics of Devon Silt under varying 
conditions of applied load and freezing rate. In each 
experiment water flow into or out of the sample, vertical 
displacement and temperature were observed and _ plotted 
against time. The program consisted of 4 series of 
experiments. The step temperature applied to all the samples 
in a particular series was the same, only the stress level 
was varied, giving a complete description of the moisture 
migration characteristics with respect to applied load for a 
given freezing rate. The step temperatures used were -5°C, - 
225°C, -1.0°C and -13°C for the NDS series, NDSB- series, 
NDSC series and NDSD series respectively. In series NDSB and 
NDSC an initial step temperature of -10°C was necessary to 
initiate nucleation, however once the sudden heave due to 
nucleation occurred,the temperature of the boundary was 
returned to the desired test temperature and the system 


allowed to equilibrate. 


To give a complete picture of the experimental 
procedure and interpretation of results, the description of 


a typical test is given in detail (NDS-7). 
The sample was prepared as described in Chapter 3 and 
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the experiment was commenced by loading the sample. Since 
the sample was prepared as a Slurry, and the top seal on the 
piston was efficient to only 1.5 Kg/cm2. The initial load 
increments were small, about 0.5 Kg/cm2 to 1.0 Kg/fcem2 to 
prevent extrusion. The excess pore pressures generated 
during each loading phase were allowed to dissipate 
completely before the next increment was applied. In test 
NDS-7, three load increments of approximately 0.73 Kg/cm2 
were applied to give a total load of 2.2 Kg/cm2. When the 
sample was fully consolidated to the desired load, the 
initial sample temperature distribution was examined by 
monitoring the readoyt from the six thermocouples placed in 
the side of the jacket of the apparatus. The requirements 
for temperature stabiity were that all the thermocouples had 
the same temperature, and that this value remained almost 
constant, any fluctuations reflecting only the small 
variation of room temperature with time. If the fluctuations 
of sample temperature were less than 0.5°C the sample was 
considered thermally stable. The temperature of the coolant 
was noted to ensure that the Sabet temperature bath had 


attained the required fluid temperature. 


Initial readings of the burette measuring pore water 
volume change and the LVDT measuring vertical displacement 
were taken and at a convenient time, eg., on the hour, the 
experiment was started. The temperature of the sample was 
continuously monitored on the strip chart recorder, and the 


readings of vertical displacement and volume change were 
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recorded manually at ten minute intervals for the first hour 
and at longer intervals, usually twenty minutes, for the 


remainder of the experiment. 


The results of experiment NDS-7 are presented in 
gcaphical form in Fig 4.1. The experiment was allowed to run 
for such time as was necessary to ensure that equilibrium 
had been achieved. Equilibriunm was defined as the 
development of a constant slope on the moisture migration 


versus time plot. 


The graphs obtained for NDS-7 were typical of the other 
tests. Considering the pore water volume change versus time 
curve, it was noted that at first a small volume of water 
was expelled. This trend lasted only for the first hour at 
which time the overall volume change was zero. After this 
point, equilibrium was quickly developed and maintained for 
two hours. The value of moisture migration rate was obtained 
from this graph and was calculated as the gradient of the 
line once equilibrium had been developed. This value was 
then used as a point on the moisture migration versus 
applied load curve, from which a more complete picture of 


migration characteristics may be oftained. 


The primary expulsion can be interpreted as a result of 
two factors. First the application of the cold temperature 
boundary condition caused the temperature of the soil close 
to the boundary to cool very rapidly. In time the pore water 


temperature dropped below the freezing point, however 
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Fig. 4.1 Sample Experimental Results NDS 7 
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nucleation did not occur immediately, as explained in 
Chapter 2,and the soil matrix close to the boundary 
continued to cool. At this point there were no excess pore 
pressures, Since to generate these nucleaction or volume 
change due to phase change is necessary. AS nucleation 
occurred the frost front propagated rapidly through the soil 
matrix that had been supercooled, which caused a distinct 
volumetric increase at the base of the sample. The resulting 
upward force against the applied load generated a pulse of 
excess pore pressure which immediately expelled some water 
and hence the pore pressure dissipated. However due to 
nucleation, suctions were being generated at the frost front 
resulting in a change of direction of moisture flow, i.e. 


the sample started to suck in water. 


The heave versus time curve for NDS-7 had a very slight 
curvature. Heave due only to segregated ice plotted as a 
straight line due to the constant rate of water 
inflow, however the non-linear advance of the frost front, 
and the corresponding non-linear heave due to freezing of 
in-situ water induced this curvature into the total heave 
behaviour. The heave rate was taken as the slope of the 


heave versus time graph once stability had been achieved. 


The temperature data from the strip chart recorder was 
plotted as in Fig 4.1. From this data a graph of depth of 
penetration of the zero degree isotherm versus the square 


root of time was plotted and the rate of advance of the 
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frost front calculated. A slight time lag between the true 
Sample temperature and that measured by the thermocouples 
due to the insulating effects of the membrane and grease was 
inevitable. It was felt that this was a minor source of 
error, since the frost front was moving sufficiently slowly 
that a lag of one minute would have a negligible effect on 


the calculation of freezing rate. 


After termination of the experiment, the apparatus was 
disassembled and the sample removed for visual inspection. 
The presence or absence of ice was noted, and then the 
sample moisture content and void ratio for both before and 
after the experiment was calculated. The void ratio against 
effective stress curve for the sample material 1S given in 


Fig 4.2. 


.2 Results of Series NDS 


The moisture migration, heave, and temperature 
characteristics of the experiments in Series A are presented 
in Figs A.1 to A.13 (Appendix A) and are summarized in Table 


pay 


The dependence of moisture migration rate on applied 
load is shown in Fig A-1. The load at which zero flow rate 
in or out of the sample occurred was 3.1 Kg/cm?. All _ the 
data points lie within acceptable limits ,however, 
experiments NDS-8 and NDS-9 experienced problems with 


temperature control of the room and the coolant bath. In 


abit ads. aoew tod pel nat | 
esiquocjsaieds oft il boxuess dead, bas, 214 tEIO" 
: <a 


7a ; | ‘ai . re 
HOY SELLS 7p hak aastdnee od? 36 i090 339 put : ae . an 


1G Go TuOF 


i P > ‘ f ‘ 
‘ 
Tha] we 74 
- ’ ~~ « 
+ ha BV eG 
ec imeagnr ) 
» t ‘ 


od7 agin One 


ine 9 otad 


3268 ses “ ig 


| ae 
. ; c \ 
ayuteteahe* nis ~SY EO ike ay ore oe in au Jaton baie 


he fnezsiy S364) 231252 Aa £23 (ar L59QkS ea" Io pou sebaey 


gideT a4 bexiabanve sie Los () xibasaqgad) ff A 


betlags no #1 a nokTAIVin oanjaton ae soushaoget 

: : vet : 

da Molt a10s doidw ts boot war +f, “+ yim 
' ? 


a , 
ats ila 


7 gs hes sats 
a re 
’ 


fed 
7 rs 
on) 


, 


Ve | 7 
' 


yan 
art 


Bibonk ret. (i Hers elqaée 
ie 


cae hd rv € raid “a as ais Ba S . att 
ai | A ; D : q | 


aA 
; i Ls 


i 
aoe i a 


iene 6 aan etea Faits 3f33 an itd Bat ve ; 


ishorlige paivon ten: $0077 12077 aaa a 


€ er 
siti. ioe s¢e hd vow asin rn aio a eae 
ra iat 
stey peassex? To foijelaaisa am 
a 
\+a@ : 
. ; ; 
Tok AI 5% une sag 10 pot tens ies ” ay 
Liti2aiv. ws daveser elie. ont Dre hn cS 
fatog ‘¢ <i 30. banauRe | Je sa _ ~— 


AY. 


ig he ha & Terr, Aue Be = En . 
52 tidy 7 -pes fooleaS = 6 W' Shak'Ectete ot aa he d 


ei ferusiiw aidgse a4). 1oJ eynus R26 378 ay tes abe 


Hg0d tC) Os 7h 


> 
» 


7 


+h 1. 
a > . 
a. a 7 
* 


bs ” 
oe 
‘bis ko 


’ 
’ 


oo - 


7 


Sgr ee 


% Ee 


ee At 


Ratio 


Void 


1.00 


90 


0 


60 


505, 


Gi. 


4.2 


0. ihe) 
Effective Stress (kg./cm? ) 


Void Ratio v's Applied Load 


50 


10.0 


— — “4 
Seen . 
ee am a, lms genes 


:  —_ as =~} 
ee 
= ls 

= a a j 

as 
= a ie és a 


. iad 
a : > 


semk, ‘entre 3 


im of 


Results of Series NDS 


STEP TEMPERATURE =-5°C 
PORE WATER 


APP : 
EX PERIMENT rare | ore MASS 
LOAD COUPLE | VOL. CHANGE 


See ee ee m.° /sec. 


0.28 2.74x10- 3.78x10° 950x10 © 6.00 x10 ° 
0.55 390 « 3.90 " 42x10" 5.42" 
2.80 4,00 (ovens 5.42 " 
280 3.90 Lol" 442 4 
2.85 390 9.00x10 ° 2.37 
3.36 4.00 769 Lgl on 
3.80 4.07 756" 25 4 
3.60 4.10 534. 4 0.81 
4.14 555 « | o7o 
4.10 476 0 =0.22 1 
332" = O35" 
250 1 - 0.42 « 
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spite of the premature termination of these experiments a 


value of was obtained but used with caution. 


In the low stress range (0-2.5Kg/cm2) the rate of 
moisture migration is particularly sensitive to applied load 
yet above the shut-off pressure it is found to become 
asymptotic to a value approximately 0.75x10!9 m3/sec. Since 
overall volume expansion was almost completely restrained at 
these high loads, this value reflected the amount of 
moisture expelled from the soil voids by the expansion of 
in-situ water on freezing as the frost front advanced 


through the soil. 


The temperature data was plotted as depth of 
penetration against the square root of time (Fig 4.3 and Fig 


4.4). Freezing rates were described in terms of X where 


X=a fC 
Where X = depth of frost penetration 
a = a constant 


t = time 
-§ te -¢ Vg 
values of 2.74x10 m/ sec to 4.3x10 m/ sec when found. 
Predictions using the simplified Stefan solution (Nixon 
~ -f r) F 
1973) of 3.78x10° m/ sec*to 4.25x10° m/ sec™ were obtained. 
The correlation between observed and predicted was best at 
high applied loads because of the lower rates of moistue 


migration in this range, a factor not considered in the 


development of the Stefan solution. 


B tests to ensure full load transfer to the sample were 
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carried out in test NOS 6 before and after application of 
some of the load increments. B values of unity were obtained 
(Figs 4.5 and 4.6) showing that full load transfer was being 


achieved and that the Samples were completely saturated. 


The total heave rate in each experiment decreased with 
increased applied load due to the reduced rate of water 
Migration. Heave due to freezing of in situ moisture was 
approximately constant in each experiment since the rate of 


advance of the frost front was similar in all the tests. 


Visual inspection of the samples after disassembly of 
the apparatus was carried out, and it was noted that the 
occurrence of segregation ice lenses, both horizontal and 
vertical, were more evident in the low applied load samples 
(Plate 1). The thickness and spacing of the horizontal 
lenses varied along the length of the Sample. Close to the 
base plate where the velocity of the frost front was 
relatively high, the lenses were thin and Garels visible. 
The lenses became progressively thicker and farther apart 
with increasing distance from the cold, side boundary. The 
increased thickness and spacing is caused by the difference 
between the rate of advance of the frost front and the rate 
of ice lens growth, which is limited by moisture 
availability, i.e. if the frost front passed rapidly through 
a soil element as in the case of the soil close to the cold 
boundary, that element would be unable to suck in an 


appreciable amount of moisture during the interval required 
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for freezing of the element. Consequently segregated ice 
lens growth at this position was subdued. However as the 
frost front progressed through the Sample and its velocity 
reduced, a larger volume of migrated water may enter the 
same volume of element in the time necessary for freezing, 


resulting in thicker ice lenses at depth. 


Vertical lensing occurred at very low stress levels. 
The negative pore pressures induced by freezing had the 
effect of increasing the effective stress on the unfrozen 
soil. The soil experienced a volumetric contraction which 
manifested itself in the form of hexagonal tension cracks 
which then filled with water and subsequently froze, civing 
the vertical lenses, refered to by MacKay (1972) as 


reticulate ice forms. 


With increasing load, lensing was progressively more 
subdued and at applied loads of 1.5 Kg/cm2 or greater was 


barely visible at any point in the sample. 


4.3 Results of Series NDS~-B 


Experimental results are presented in graphical form in 
Fig A.-14 to A.20 (Appendix A) and summarized in Table 4.2. 
Fig A.14 summarizes the moisture migration characteristics 
due to a change in applied load. In test NDSL-T-1 the sample 
was consolidated to 5.0 Kg/cm? and freezing was permitted to 
continue for a much longer time period than the other tests, 


and an important characteristic was noticed (Fig A.19 and 
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Results of Series NDS-B 


STEP TEMPERATURE =—25°°C 


APPLIED o< THERMO - HEAVE PORE WATER 
EXPERIMENT 
LOAD COUPLE RATE VOL. CHANGE 
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vse. [| m/vsec | m/sec. | m/sec. 


NDS B-I 70x10” POOL 8.33x10 ° 6.067x 10° 
NDS B-2 mona 2.80 350 1 
NDS B-3 290 


NDS B-4 : 296 


NDS L-T-| ————— 300 


Table 4.2 
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Fig A.20). Expulsion of pore water was the Orimary moce of 
behayiour, however after 550 minutes expulsion stopped and 
water was attracted into the sample. This trend continued 
for approximately 1400 minutes ( 1 day) and the test was 
terminated after this time. This change in moisture transfer 
characteristics shows very clearly the limitations of ee 
assumption that short term behaviour accurately represents 
the long term behaviour of a freezing soil. The long tern 
behaviour shows that the shut-off pressure is greater than 
5.0 kg/cm? however using short-term data a shut off pressure 
of 4.0 Kg/cm? is obtained. It is felt that the long tern 
data exhibits the true mode of behaviour and that the 
moisture migration Ay See eee obtained in the short 


term tests must be interpreted with care. 


On visual examination of the samples, ice lensing as 
described in Chapter 4.1 was noted, however the frost front 
had not penetrated as far into the samples as in the NDS 
series and consequently only the fine closely spaced lenses 
ty pical of those formed close to the cold side boundary were 


found. 


The magnitude of total heave rate decreased with an 
increase in applied load reflecting the decreasing rate of 


moisture migration at higher stresses. 


The difference between values predicted by the Stefan 
solution and observed values became greater at this step 


temperature due to the errors involved in neglecting the 
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effects of migrated moisture and the error in assuming the 


initial sample temperature to he 0°c. 


4.4 Results 


The results of the experiments in Series NDS-C are 
presented an graphical form in Figs As21 "to  A.23 (Appendix 
A) and summarized in Table 4.3. The shut-off pressure was 
found to be 4.2 Kg/cm2. The same limitations due to the 


effects of short term testing apply here as in Series NDS-B. 


To check the validity of the use of short duration 
tests for the low applied load samples, tests NDS C-1 NDS C-2 
and NDS C-3 were permitted to run for 640 minutes, 640 
minutes and 1420 minutes respectively, to ensure that the 
behaviour in the long term was consistent with that assumed 
from the short term tests. The long term data points fell on 
the projected lines (Figs A.28 A.29 and A.30) confirming 


this assumption. 


The magnitude of the rates of migration of water are 
lower for a given applied load than in the _ previous 
series,the implications of which will be discussed in detail 


later. 


In this series of experiments the frost front became 
stationary thus making the calculation of values 
meaningless and since the frost front did not pass any of 


the thermocouples derivation of rate of advance by this 
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Results of Series NDS-C 


Ser ie MPER ATURE = 15.6 
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Table 4.3 
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method was impossible. 


The total heave rates were also reduced in comparison 
with the previous series due to ‘the decreased rates of 


moisture migration. 


Lensing was noted in the samples at low applied load 
es pecially NDS C-1 which sucked in 6.0x10* m3 to a 
Stationary frost front, building up a comparatively thick 


layer of ice. 


The results of the experiments in Series NDS-D are 
presented in graphical form in Figs A.29 to A.32 (Appendix 


A) and summarized in Table 4.4, 


The shut-off pressure was found to be 4.0 Kg/cm2 (Fig 
A.29) The graph of rate of pore water volume change against 
applied load has the same shape as the curves obtained in 


the NDSA series and the NDS-B series. 


/The volume change against time curve in experiment NDS 
D-1 exhibited continuous curvature due to the very rapid 
advance of the frost front. Thus to obtain a value of BO ete 
the gradient of the curve was taken at a time when the frost 


front was well inside the middle third of the sample. 
Lensing was not visible in any of these samples. 


The total heave rates were similar to those in the NDS 
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series, in response to the similar rates of moisture 


migration. 


Correlation between predicted and observed values is 
not good, ranging from 24% error to 30%. It is thought that 
the extreme value of step temperature applied may have 
caused problems such as heat conduction into the sample 
through the walls due to the high temperature gradient and 


induced these errors into the systen. 


4.6 Conclusions From Observations 


In all of the four series of experiments, an increase 
in applied load was accompanied by a decrease in rate of 
moisture migration due to a combination of two factors. 

1) A decrease in the total potential gradient of the pore 
water with increase in applied load, 
2) A decrease in the soil permeability due to increase in 


applied load. 


At high applied loads the moisture migration 
Characteristics obtained from short duration tests did not 
not reflect the long term behaviour, as found in test NDS L- 
T-1. Consequently the values of shut-off pressure obtained 
in each series must be considered inaccurate. Hence the 
evidence obtained in this test program which showed that the 
Shut-off pressure was independent of temperature 
distribution is inconclusive since it is felt that the true 


value of shut-off pressure was greater than 5 Kg/cm?, and 
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not about 4.0 Kg/cm2 as indicated by the short term tests. 


Comparison of the moisture migration characteristics 
obtained at low applied loads in short duration tests is 
valid, since long tern testing in series C showed that the 
behaviour of freezing soil at low applied load remained 
unchanged with time. Comparison of the moisture migration 
rate against applied load curves for series NDS, NDS B, and 
NDS D show the initial portion (0 Kg/cm2 to 3.0 Kg/cm2) to 
be identical, however the series NDS C curve falls well 
below these. Considering an experiment from each series at a 
given stress level, the permeability is the same in each 
case, hence the reduction in water flow in to the sample 
must be a ene of a reduction in the suction generated at 
the frost front. The cause of this reduciton wili:be 
discussed in Chapter 5. These results imply that the suction 
generated at very low heat energy flow is small, increasing 
as the heat flux is increased, until a plateau value of 
suction, independent of heat flux, is maintained. This value 
Of Psiiction is probably a function of soli particle size 


distribution and stress history. 


The magnitude of the total heave ratio obtained by 
measurements in each experiment and the segregational heave 
rate obtained from the moisture migration characteristics 
were plotted against applied load. (Figs 4.7 to 4235) 2 itis 
clear that both the total heave rate and the segregational 


heave rate were reduced by an increase in applied load. The 
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ratio of these two components of heave was also plotted 
against applied load. The effects of assuming the.expulsion 
behaviour at high applied load also renders the upper 
portion of the curve in doubt, however it is clear that at 
low applied loads the ratio of segregational heave to total 
heave decreases with an increase in applied load. At the 
higher step temperatures, the ratio of the heave components 
1s approaching unity which implies that fully efficient ice 
HONS @eng 1s Occurring; Lf “this ratio is, defined “as “the 
efficiency of the lensing system, the graph of efficiency 
against the square root of applied step temperature (Fig 
4.16) shows that the efficiency of the system decreases as 
the advance of the frost front through the sample becomes 


more rapid. 


The major limitation of the experimental data is the 
validity of the high applied load results. Experiment NDS L- 
T-1 demonstrated the change in behaviour between short 
duration and long duration testing at high stress. 
Consequently in any further research work the importance of 


the duration of the test must not be underestimated. 


In the anaylsis of the test results values of suction 
developed at the frost front were derived indirectly from 
the value of moisture migration rate observed in 


experiments. The main limitation of this derivation is that 
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the value of suction calculated is very sensitive to the 
permeability of the ‘soil. ‘The determination of the 
permeability of the soil is in itself an inexact procedure 
since permeability is very sensitive to the set-up procedure 
and aandling of the apparatus. The value of permeability 
used in the calculation of the suctions developed was 
obtained from a constant head permeability test using the 
Permode apparatus thus minimizing any inherent errors as 


much as possible. 


Direct measurements of pore water pressure during 
freezing were not taken since it is necessary to have a 
closed system drainage condition for the measurement of 
water pressures at the end of the sample. Since the 
ex perimental program was investigating open system 
behaviour, and it is felt that the application of a closed 
system boundary condition alters the freezing conditions, 


these measurements were not attempted. 


La on 7 we vv 


eit oF sVistenes YISV 4 


: . 7 
ant +0 nortantmietab aa 


anubeaorg 2 2Aisn i as [Lec v9#) dé, et bioew oat 
sasyheasoaq qu-ya2 ed? 67 evi teas? yTev at AL 
sd? set Se7T ERTS ott bail enibe 


7 
snofssuc Ont 36 cceeenboten | ede a 


gf#itidesmysg. ae Buse 


“ G chet fou nt 
26V DSUcCl.ov>oY : 


(2a si{idsomteq' btied tasetanoo & word t 


ak =in2ae f#0ersdal ys ehisieigre |) eont auseiegys. 


it 
-Stdinaod 68 
7 
= 4 os aus : “mm NoAG 4 Te 
pocztl Swvessiq 797 6v TOG 70. 220SGe zen sy a u ta 
ingigla lh oe 
& 9vV5A O28 YAbeaqaVSn es rs sonia as¢g@set 2on e798 parse x 


7 
. * & Z — 
Jo #0ens ye san ay’ 107 AoOiPihaos habia asteys i. 


esate ,.2iaane 4 Fe t ic ent Ph agtiueewa( 


qa tvave nega Paste = sit ce wstpoiwd Sen 

j ‘an we nae 

saols io Yottedilage eft ted? Gia al nae’ +18 
} a 


of398 13 oai4 2iad ls noobs Loa 
: eA, 


ae ee sort ate) 


ie fl j 
fafaget*s to ease etre TERR Se 


When subzero temperatures are applied to the surface of 
a soil, freezing will be induced. If the surface temperature 
history is known then a solution for the freezing rate and 
temperature distribution in the soil with time may be 


calculated. 


If the surface temperature is held constant at some 
value below 0°C, then the Neumann solution (Carslaw and 
Jaegar, 1947) is capable of describing the rate of freezing 
by the equation 

X = ave See 


Where ag 


depth of frost penetration 

a& = constant, depending on the thermal 
properties of the soil, and the ground and 
surface temperatures 


t = time. 


If the ground surface temperature varies with time, the 
preceding method may no longer be capable of modelling 
accurately the soil behaviour. For this reason a numerical 


solution must be adopted to determine the rate of freezing. 


The governing differential aquation of one dimensional 


conductive heat transfer from Ho, Harr and Leonards (1970) 
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Where 6 = temperature 
C (68) = volumetric heat capacity of the 


soil/ice/water mixture 


K(®) = thermal conductivity of the soil 
L = latent heat of the soil 
W = water content 


Ya= dry density 


W,(@) = unfrozen water content. 


Writing the above in finite difference form assuming a 


constant discrete interval in the x-direction gives 
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Where i and j are the finite difference subscripts in the x 
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those 


must be satisfied, and this limitation on the time step t 


is the chief drawback to this numerical method. 


The input data necessary is KurkXeePs,O,He Ts, Tyo Wi 


Where K,,K, = thermal conductivities of the fully thawed 


and fully frozen soil respectively. 
ey?) define the unfrozen 


content/temperature relationship 


equation 


Wy = (Ptexp (Q0+R)) 7100 SP 
Where R = in(100-P) 
P = represents the terminal 


water 


by the 


water 


content at a very low temperature as a 


percentage of the total water content. 


Q = defines the curvature 
relationship. 
H = height of the soil considered 


Ts, and T, are the surface and ground 


temperatures respectively 


W is water content expressed as a 


percentage. 


The units used are centimeters, grams, seconds, 


Wu (8) 


degrees 
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Centigrade. 


A listing of the program and a sample output is given 


in Appendix B. 


The results of.depth of frost penetration against the 
square root of time from the computer program are presented 
ieee) 9.) ands Fig@ 5.22 —These | figures represent the 
propagation of the O°C isotherm for various configurations 
of Step temperature and applied load. values were 
calculated for each of these and compared to the e values 
obtained in the corresponding experiment. Figs 5.3 and 5.4 
show the comparison of the values predicted with the 
values observed for applied loads of 0.55 Kg/cm2 and 3.0 
Kg/fcm2 respectively. Comparison of the predicted and 
observed is good however limitations to this development 
exist. As the frost front propagates through the sample and 
heat is extracted from the unfrozen soil a finite quantity 
of heat is conducted through the top piston providing an 
extra source of heat energy. Hence as the frost front 
advances through the soil the end effects become 


increasingly pronounced. 


In the development of the governing differential 
equation the flow of migrated moisture was neglected since 
there was no composite method of predicting water flow to 
the frost front and its effects in the heat conduction 
equation. At low stresses where the volumetric flow rate of 


water to the frost front may be of considerable magnitude 
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Frost Penetration (m.x 10-2) 


Depth of 


=2 
Penetration (m.x10  ) 


Frost 


Depth of 


2.0 


3.0 
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10 
VTime (mins.'“7) 


10 
Time 


( mins {/# 
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& = 0.555 kg./em? 


K-19 = 1.485 x 1075 msec. 
o -25= 217 x 10% m./ eee. 
XK 597293 «1075 m./vsec. 
A 13.07 4-42 x 10-5 m./ Vsec. 


O= 30kg Jem? 
Oven 5 =herel x107°m./ VSeC. 
SL 59 =2.95 x10 5m./ sec. 
OX 1302 4.461075 m./ V5ee 
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inaccuracies will be inherent in the calculations. From the 
moisture migration data observed in the test program at way 
be noted that at lower heat flux removal rates the heat 
energy brought to the frost front in the form of latent heat 
of fusion may be sufficient to be of consideration in the 


heat balance equation. This is another limitation placed on 


the use of the progran. 


A combination of these factors, viz., low applied load 
and a small negative step temperature,. may induce 
considerable inaccuracies and values of depth of frost 


penetration should be used with caution. 
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The data from the experimental program demonstrated 
that a potential gradient was induced through the unfrozen 
zone during freezing causing water to flow into the sample. 
The moisture migration rate was found to be a function of 
both applied load and heat flux removal. The affect of 
increasing the applied load was to decrease the rate of 
moisture migration. The behaviour may be interpreted in 
terms of two factors: 1) decreased permeability of the soil 
and 2) a decrease in potential gradient through the soil. As 
the applied load is increased the permeability of the soil 
decreases reflecting the decrease in void ratio. Thus even 
if the potential gradient through the soil was independent 


of applied load, the moisture migration rate would still 
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decrease with applied load due to the effects of soil 
permeability alone. However, the potential gradient also 
varies cue to a change in applied load as predicted by the 
Kelvin equation and demonstrated by Penner (1959). The 
interaction of both these factors makes investigation of the 
actual value of suction developed due to a change in applied 
load impossible, without recourse to direct measurement of 


suctions in experiments using special techniques. 


The decrease in rate of moisture migration due to a 
decrease in heat removal at low step temperatures is not 
affected by changes in permeability since the water flow 
rates are compared at the same applied load. Consequently 
the decrease in flow rate must be solely a reflecton of a 


decrease in potential gradient. 


From the moisture migration characteristics, values of 
the suction generated at the frost front may be calculated 
if the permeability of the soil is known. The soil structure 
is assumed to be incompressible, hence no moisture transfer 
due to consolidation will occur. The development of a 
constant value of suction at the frost front will induce a 
linear pore pressure distribution through the soil sample 


ie. the hydraulic gradient is a constant. 


Using D‘tarcys law, 
Ves Ka Sit 


velocity of the water 
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i= permeability 
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1 = hydraulic gradient 


The hydraulic gradient may be back calculated and hence 
the suction generated at the frost front. 
= KiA Sao 


Where area of sample 


rate of moisture migration 


PIE > lp. 
it 


Calculation of the magnitudes of suction developed are 
given in Table 5.1. These results indicate that at low heat 
flux the suction developed is a function of heat energy 
flow, however after a certain level this dependency stops 
and the suction reaches a maximum plateau value which cannot 
be exceeded. The factors limiting the magnitude of the 
Suction are probably pore size (itself a function of the 
grain size distribution) and cavitation effects in the pore 


water. Fig 5.5 shows this behaviour in graphical form. 


This is a mode of behaviour which may well embrace the 
findings of two previous studies. Beskow (1935) demonstrated 
that the heaving rate was not influenced by the rate of 
frost penetration however his atee ee temperatures were 
in -the range of -3°C to -10°C, confirming the ‘plateau’ 
concept deduced from this investigation. Penner (1959) 
showed that the heave rate was dependent on heat flux. It 
was demonstrated that increasing the rate of heat flow 
induced higher moisture migration rates and consequently 


higher heave rates. It is felt that the marriage of these 
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two modes of behaviour may give the complete picture of the 


effécts of heat flux. 


The physical interpretation of the variance of the 
suction developed may be made with respect to the concept of 
an energy budget. At low step temperatures where the rate of 
heat flow is low there is a limit to the amount of water 
which can be attracted to the frost front. This is due to 
the fact that the water flowing towards the front carries a 
certain amount of heat energy in the form of latent heat of 
fusion. This heat flux is in addition to that brought to the 
frost front by conduction and this total cannot exceed that 
removed from the frost front by conduction through the 
frozen zone. Thus if the magnitude of heat flux removal was 
very low a reduction in mnoisture migration rate would be 
anticipated, and was in fact observed in the experimental 


program. 


The validity of the linear pore pressure distribution 
assumption is of interest. Since the soil structure is not 
an incompressible medium, the pressure gradient through the 
sample will not be constant until a finite amount of 
consolidation has occurred. Also the dynamics of the frost 
front further complicate the problem. The effect of these 


factors is to introduce curvature into the suction gradient 


{Fig 5.6). 


This curvature is most pronounced immediately after 


nucleation. The development of equilibrium conditions is 
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Fig. 5.6 Schematic Diagram of Pore Pressure 
Distribution of the Soil Sample 
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also hindered by the advancing frost front. As the frost 
front moves through the soil consolidation has to occur 
through the sample, and equilibrium will only be developed 
if the frost front remains stationary sufficiently long for 


dissipation of excess pore pressures to occur. 


The application of a very low step temperature would 
cause a rapid advance of the frost front and hence would not 
permit the development of equilibriun conditions. 
Consequently a high suction gradient would exist at or near 
the frost front lending support to the observation of MacKay 
(1974) that during rapid freezing of freshly exposed fine- 
grained sediments in the high Arctic, water is drawn from an 


area close to the freezing front. 
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CHAPTER © 


a) There is a relationship between the rate of moisture 
migration into a soil on freezing and the applied load. The 
rate of moisture migration falls sharply with an increase in 
applied load until the situation of no flow into or out of 
the sample occurs. The values of shut-off pressure obtained 
from the experimental program cannot be compared with each 
other since the high applied load data (greater than 3.5 
Kg/cm2) is inconclusive due £0 the effects of short term 
testing. 

b) To obtain the true moisture migration characteristics at 
high applied loads, long duration testing must be employed 
Since it was shown that migration behaviour may change after 
some time at high loads. This is thought to be the due to 
the effects of secondary consolidation. 

c) The rate of moisture migration is a function of the rate 
of heat energy removal. As the rate of heat energy removal 
increases the rate of moisture migration also increases 
until a maximum ‘plateau' vaiue is reached. The suction 
causing this maximum value of rate of water flow is probably 
a function of grain size distribution and stress history. 

d) The temperature distribution during freezing may be 
calculated using one dimensional heat conduction equations. 


A finite difference formulation of the equations was found 
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to predict toa high degree of accuracy the advance of the 


0°C isotherm through the soil samples. 


ie 
ee) 


Equipment Performance and Procedure 


Sa 


The apparatus, although not designed for work with 
freezing soils, performed well. The problems encountered 
with other freezing cells which employ different sealing 
systems were avoided by the use of the O-ring and groove 
construction. Leakage did not occur since the , load 
increments applied during the consolidation phase were 
small. Full saturation of the soil samples was obtained 
using the vacuum and vibration technique,however care must 
be taken when placing the sample in the apparatus to avoid 
trapping pockets of air between the soil and the side walls 
of the cylinder. Friction was also virtually eliminated and 
full load transfer to the samples was obtained. B tests were 
carried out on various experiments to check for saturation 
and load transfer. The data obtained from the experimental 
program indicates that the overall performance of the 


apparatus was good. 


Two major improvements could be made to the apparatus. 
The use of the thermistors set into the walls of the 
cylinder instead of thermocouples would give much better 
temperature data. At times the data obtained from the 
existing temperature measurement system appeared erratic and 


unstable, and interpolation between points where the thermal 
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~ 


data appeared stable had to be carried out on occasions. ~ it 
is tert that the increased stability and accuracy 
anticipated by the use of thermistors would be a major 
improvement. Another problem was that of supporting the 
slurry while the cylinder was removed during the 
installation of the piston. This problem may be eliminated 
bypcuttang ‘the cylinder) |) jackets’ at’ >the’ required” indtial 
sample height. With this configuration the bottom portion of 
the cylinder may be used to support the slurry during 
placement. If the sample height corresponds to the top of 
the cylinder the piston may be installed without the 
necessity of removing the cylinder jacket. Having installed 
the piston, the upper portion of the cylinder may be placed 


and assembly completed. 


The test procedure as set out in Chapter 3 gave _ good 
results showing the trends of behaviour due to the varying 
conditions of applied load and heat flux removal. There may 
however be some anomaly between the results obtained and the 
actual moisture migration characteristics at high applied 
loads. Although the apparatus functioned adequately it is 
felt that the interpretation of the results may have been 
incorrect in some cases. In all the experiments the duration 
of testing depended on the shape of the moisture migration 
characteristics of the experiment Bigg = only after 
equilibrium had been developed and maintained for some time 
was the test terminated. In the case of the high applied 


load tests where the mode of behaviour was considered to be 
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91 
expulsion some difficulty may have arisen. In these cases 
the moisture migration characteristics which were assumed a 
have been in equilibrium may have in fact been a quasi- 
stable situation, i.e. in the short term the results 
indicated that a constant rate of expulsion had _ been 
achieved and the test was terminated, yet perhaps if the 
experiment had been permitted to run for a considerably 
longer period the mode of behaviour may have changed and 
water may have migrated into the sample. Thus the initial 
expulsion found in these tests may have exhibited such 
slight curvature that in the short term the moisture 


Migration characteristics plotted as a straight line. 


Another factor which may have had some affect on the 
high load tests was that of secondary consolidation. In each 
test complete dissipation of excess pore pressures oe 
ensured by plotting the vertical compression of the scil 
sample against the logarithm of time. If the curve indicated 
that the primary consolidation was completed, freezing of 
the sample could then be commenced. Hence when freezing was 
first started the sample would still have been undergoing 
secondary consolidation. In most engineering situations the 
magnitude of secondary consolidation 1s negligible in 
comparison with that of primary consolidation. However in 
this case where the rate of pore water volume change was 
very small, the effects of secondary consolidation may have 
been comparatively important foryvtheeficstesarnougs sos Etis 


test, during which time the sample was trying to,.establish 
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~ 


equilibrium conditions. Hence the very small volume of water 
expelled due to secondary consolidation may have also 
contributed to a misinterpretation of the short term data. 
Consequently the data points at high applied loads which 
were plotted as pore water expulsion RayRnotabe feorrect.. 17 
is felt that if in the long term the sample did attract 
water, the magnitude of the rate of pore water volume change 
would be very small and the trend of decreasing moisture 


transfer with increased applied load would be continued. 


This is a point of interest in further research, 
however the author is confident that all the other 
ex periments, in which the mode of behaviour was attraction, 
had developed true equilibrium and that short term behaviour 
can represent accurately the long term freezing behaviour of 


soils. 


The problem of nisinterpretation of the expulsion data 
does not affect the finding of lower moisture transfer rates 
for a given applied load due to a decrease in heat flux 


removal. 


6.3 Recommendations 


The distribution of pore pressure through the sample is 
of major importance in developing a comprehensive solution 
to the frost heave problem. The change in this distribution 
with time due to an advancing frost front, and the time 


taken for development of equilibrium conditions for a 
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93 
stationary frost front are of much interest. It is felt that 
any deformation of the unfrozen soil and the change in pore 
pressure distribution with time are functions of the 
conventional soil mechanics parameters ay, the*co=eiErcient 
of volume compressibility and Cy, the co-efficient of 
consolidation. If the magnitude of the suction generated at 
the frost front is known for a given set of conditions of 
applied load and heat flux removal, some method of 
predicting the pore pressure distribution with time may he 


developed. 


A finite difference approach may be of heip in 
predicting the pore pressure distribution in the non- steady 
state condition of an advancing frost front. Calculation of 
the pore pressure distribution through the sample with each 
time step would give a theoretical prediction which could be 


examined in the light of experimental findings. 
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Step Temperature -13° C 
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APPENDIX B 


PROGRAM LISTING AND SAMPLE OUTPUT 


Ne PE 


OOWYAMN SE WH oo 


83 


91 FORMAT (® 


42 FORMAT (*'O®,*THE TINE STEP 1S*,F10.3,! 


REAL U (200) , UNEXT (200) 
COMMON P,Q,R8,TCT,TCF,W,GAND 


t READ (5,2) TCT,TCF,P,Q,H,TS,T1,W 
2 FORMAT (8F10.0) 


WRITE (6,83) FCT, TCP,H 


134 


PORMAT(*1°,//* THERMAL CONDUCTIVITY IN THAWED SOIL IS*,F10.5/ 
** THERMAL CONDUCTIVITY IN THE FROZEN SOIL IS,°,P10.5,* CAL/SEC-CK- 
*DEGeC.°/*° THE HEIGHT OF SOIL COWSIDERED IS,*,Ff8.2,°% CH*) 


WRITE (6,91) TS,T1I,@ 


",°THE CONSTANT SURFACE TEMP IS*,F(0.5/°THE INITIAL 


*GROUND TEMPERATURE IS*,P10.5/*THE WATER CONTENT IS°,F1%0.5) 


R=ALOG (100.-P) 

CALL PLOTUF (P,Q,8) 

I1=0 

N=20 

DX=H/N 
GAMD=2.7/ (1. ¢8%2.7) 
CCF=TCF 

IP(TCT.GT.TCF) CCF=TCT 
DT=GAMD*DX#DX# (0. 2+0. 548) / (2. 2%CCP} 
DTO=0.2«DT 

DTF=DT 

WRITE (6,42) DT, DX 


*/) 

DO 3 I=1,N 

U (I) =TI 
TIME=0. 
TOL=0.02 
U(N+1)=U (N-1) 
BET=DT/ (DX«DX) 
mo) Gee) 0 

T=U (1) 

TP=0 (I+ 1) 

IF (I. EQ. 1) GO TO 6 
TM=U (I-1) 

GO TO 7 

TM=TS 


SEC. 


A=2.4#BET/ (CO(T) «(4e/TC(TP) + 16/TC(T))) 


IF(I.GT.1) GO TO 44 
B=2. *BET*TC (TS) /CO (T) 
GO TO 45 


B=2.«BET/(CO(T) «x (i-/TC(T) + 1./TC(TM))) 


AB=AtB 
IF(AB.-LT.1-) GO TO 46 
WRITE (6,47) 1 


FORMAT (*0®,* THE*,I3,'TH NODE IS UNSTABLE, 


1 BE HALVED...°) 

DT=DT/2.0 

GO TO 4 

E=79.6%W*GAMD/CO (T) 

D=AeTP+ (1.-A-B) «T+ Be TMtEx WU (T) 


KOUNT=0 

SOL=T 

STORE=SOL 
P=SOL-D+E«WU (SOL) 

IF (SOL.LT.0.) GO TO 15 
IF (T.LT.0.) GO TO 50 


COMMENCE THE ITERATION TO PIND U(I,J+1)-~- 
LET THE INITIAL APPROXIMATION BE U(1,J) oon 


DX IS*,F8.3,° CM.°*// 


AND THE TIME STEP WILL 
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50 


15 
16 


12 


13 


70 
18 


Eo 


FDASH=1. 
GO TO 16 

PDASH=1. ¢ Ex (WU (SOL) -W0(T)) /(SOL-T} 

GO TO 16 

PDASH=1.+EeQ EXP (Q*SOL¢R) /100. 
SOL=SOL-F/FDASH 

DIF F=ABS (SOL-STOBRE) 

KOUNT=KOUNT+ 1 

IF (KOUNT.LT.20) GO TO 9 
SOL=0. 5% (SOL+STORE) 

GO TO 11 

IF(DIPF.GT.TOL) GO TO 8 

UNEXT (I) =SOL 

CONTINUE 

TIME=TIME*D? 

HOURS=TIME/3600. 

SQHR=SQRT (HOURS) 

DO 12 I=1,N 

U (1) =UNEXT (I) 

IF(II-LT.100) GO TO 70 

II=II-100 

WRITE (6,13) TEME, HOURS, SQHR, (U (I) »I=1,N) 
FORMAT (* ©,*TINE=*,F12.1,° HOURS=",?10.2,° ROOT TIMNE=*,P10.3/ 
1(10F10.4)) 

GO TO 18 

TI=II+1 

IF(U(10).LT.0) stop 

GO TO 4 

END 

REAL FUNCTION TC(T) 


SSeSo THIS FUNCTION EVALUATES THE CONDUCTIVITY AS A FUNCTION OF 
aoe TEMPERATURE..-.~ 


HERE EREKEEKERSE KEE KERR SSERER EEK EERE SHE 


COMMON P,Q,8,TCT,TCF, W,GAHD 
IF (T.LT.0.) GO TO 1 

WU=1. 

GO TO 2 

WU= (Pt EXP (Qe T*+R)) /100. 
TC=TCF+HU« (TCT-TCF) 

RETURN 

END 

REAL FUNCTION CO(T) 


-THIS FUNCTION EVALUATES THE APPARENT VOLUMETRIC SPECIFIC HEAT 


-OF THE SOIL... 
~KEKEREKEEREKEEKEEESE EERE EREEE KEE SESE EERES SEE 


COMMON P,Q,8,TCT,TCF,W,GAMD 
IF (T.LT.0.) GO TO 1 

wU=1. 

GO TO 2 

WU= (P#EXP (QuT+ R)) 7100. 
CO=GAMD» (0. 2+0. 5 Wu (1.+WU) ) 
RETURN 

END 

REAL FUNCTION WU(T) 


-THIS FUNCTION BVALUATES THE UNFROZEN MOISTURE CONTENT OF THE SOIL 


-AT THE TEMPERATURE °°T*%® weeeee 
—HEEEEEK EERE ERE EERE EERE EEE EET ETE OE EE ERT ERE HE 


COMMON P,Q,R,TCT,TCF,W,GAMD 
IP (T.LT.0-) GO TO 1 
WU=1. 
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126 GO TO 2 

127 1 WU= (P+EXP (Qe TeR)) 7100. 

128 2 RETORBN 

129 END 

160 SUBROUTINE PLOTUF (P,Q,R) 

161 Se THE SUBROUTINE PLOTS THE UNFROZEN MOISTURE CONTENT AGAINST TEMPERATURE 
162 INTEGER LINE(100) , AST,DOT,BL 

163 BL= 1077952576 

164 AST=15477 14624 

165 DOT= 1262501952 

166 WRITE (6,1) 

167 1 FORMAT (°0*,*PLOT OF THE ONFROZEN MOISTURE CONTENT AGAINST TEMP. ¢/ 
168 1* TEMSP.*) 

169 DO 2 I=1,60 

170 2 LINE (1)=DoT 

171 WRITE (6,3) (LINE(L) ,1=1,60) 

4172 3 FORMAT (* #,9K,°0%,24X,°50%,23X,°100",10X,*PERCENT M.C.°/60A1) 
173 TEMP=-0.25 

174 6 UFAHC=PtEXP (Q«xTEMPtR) 

175 LL=UFMC/2.*10. 

176 DO 4 IT=1,60 

177 4 LINE (1) =BL 

178 LINE (10) =DOT 

179 LINE (LL) =AST 

180 WRITE (6,5) TEMP, (LINE(L) ,1=10,60) ,UFAC 

181 5 FORMAT (* °,F6.1,3X,51A1, 10X,F10.2) 

182 TENP=TENP-0.25 

183 IF (TEMP.GE.-5.00) GO TO 6 

184 RETURN 

185 END 


END OF FILE 
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END OF FILE 


THERMAL CONDUCTIVITY IN THAWED SOIL Is 


THERMAL CONDUCTIVITY IN THE FROZEN SOIL IS, 
THE HEIGHT OF SOIL CONSIDERED IS, 


THE CONSTANT SURPACE TEMP IS 


THE INITIAL 


THE WATER CONTENT IS 


GROUND TERPER 
35.53000 


6.00 
-5.00000 
ATORE IS 3 


0.00289 
Ca 


«44000 


OPLOT OF THE UNFROZEN MOISTURE CONTENT AGALNST TERP. 


TEMP. 
0 


50 


Se OS SSO RSS MSOC ECC ECS SSA CUU CEC CESS We CCC EEECEVecceeesobe 


-0.3 . 
-0.5 . 
-0.8 : 
-1.0 . 
-1.3 . 
SoS) . 
-1.8 - 
-2.0 . 
iS) . 
=2e2 . 
-2.8 - 
-3.0 . 
-3.3 : 
-3.5 : 
-3.8 - 
-4.0 . 
-4.3 ° 
“4.5 - 
-.8 : 
-5.0 : 
OTHE TIME S 


TIME= 
-3.6559 
3.2044 
TIME= 
-4.0925 
2.6294 
TIME= 
-4.2391 
2.1557 
TINE= 
-4.3466 
1.7812 
TINE= 
-4.4140 
1.4717 
TINE= 
-4.4657 
1.1955 
TIME= 
-4.5059 
0.9681 
TIME= 
-4.5381 
0.7593 


TINE= 
-4.5651 
0.5918 
TINE= 
-4.5881 
0.4438 
TINE= 
-4.6078 
0.3236 
TIME= 
-4.6252 
0.2269 
TINE= 
-4.6406 
0.1468 
TIAE= 
-4.6544 
0.0901 


TEP IS 4.127 SE 
416.9 HOURS= 
-0.7951  -0.0078 

3.2926 3.3505 
833.7 HOURS= 
-2.2630 -0.4119 

2.8143 2.9640 

1250.5 HOURS= 
-2.7115 -1.2006 

2.3699 2.5545 
1667.4 HOURS= 
-3.0360 -1.7300 

1.9980 2.1884 

2064.2 HOURS= 
-3.2391 -2.0676 
1.6754 1.8573 
2501.% HOURS= 
-3.3949 -2.3265 
1. 3880 1.5599 
2917.9 HOURS= 
-3.5158 -2.5273 
1.1397 1.2945 
3334.8 HOURS= 
-3.6128  -2.6886 
0.9178 1.0587 
3751.6 HOURS= 
-3.6939 -2.8237 
0.7278 0.8514 
4168.4 HOURS= 
-3.7631 -2.9387 
0.5651 0.6745 
4585.1 HOUBS= 
-3.8224 -3.0376 
0.4257 0.5193 
5003.7 HOURS= 
-3.8747 -3.1245 
0.3147 0.3945 
5418.3 HOURS= 
-3.9210 -3.2016 
0.2189 0.2855 
5834.9 HOURS= 
-3.9623 -3.2705 
0.1896 0.2041 


* 
* 
* 
* 
® 
* 
Ce. DK IS 


0.12 ROOT 
0.7382 
3.3872 

0.23 ROOT 
0.1097 
3.0828 

0.35 ROOT 

-0.1249 
2.7104 

0.46 ROOT 

-0.5079 
2.3574 

0.58 hoot 

-0.9392 
2.0164 

0.69 ROOT 

-1.2826 
1.7103 

0.81 ROOT 

-1.5539 
124312 

0.93 ROOT 

-1.7747 
1.1831 


1.04 ROOT 
-1.9613 
0.9612 
1.16 ROOT 
-2.1202 
0.7712 
1.27 ROOT 
-2.2573 
0.6029 
1.39 ROOT 
-2.37862 
0.4656 
1.51 ROOT 
- 224855 
2.3456 
1.62 ROOT 
2.5814 
0.2530 


* 


0.300 CK. 


TIME= 
1.3690 
3.9096 

TIAE= 
0.6135 
3.1746 

TIME= 
0.2671 
2.8388 

TIME= 

-0.0092 
224904 

TIME= 

-0. 1356 
2.1521 

TIME= 

-0.3852 
1.8385 

TINE= 

-0.6599 
1.5487 

TISE= 

-0.9105 
1.2898 


TIAE= 
- 1.1330 
1.0561 
TIME= 
-1.3257 
0.8545 
TIME= 
-1.4951 
0.6756 
TIAE= 
-1.6460 
0.5271 
TINE= 
-1.7804 
0.3981 
TIME= 
-1.9013 
0.2955 


* 


0.340 
1.8912 
3.4228 

0.481 
1.0670 
3.2436 

0.589 
0.6469 
229413 

0.681 
0.3314 
2.6027 

0.761 
0. 1668 
2.2637 

0.834 

-0.0090 
1.9441 

0.900 

-0.0902 
1.646% 

0.962 

-0.2469 
1.3780 


1.021 
-0.4317 
1.1351 
1.076 
-0.6153 
0.9237 
1.129 
-0.7933 
0.7363 
1.179 
-0.9585 
0.5784 
1.227 
1. 1094 
0.4421 
1.273 
-1. 2480 
0. 3313 


100 


0.0049% CAL/SEC-CM-DEG.C. 


2.3123 
3.4302 


124717 
3.2933 


1.0032 
3.0191 


0.6621 
2.6896 


0.4619 
2.3508 


0.2537 
2-0267 


0.1409 
1.7227 


-0.0016 
1.4472 


-0.0459 
1.1973 


-0.1409 
0.9782 


-0.2658 
0.7843 


-0.3966 
0.6190 


-0.5321 
0.4770 


-0.6666 
0.3596 


PERCENT #.C. 


86.73 
76.39 
68. 34 
62.07 
57.19 
53.39 
50.43 
48.12 
46.32 
44.93 
43.84 
42.99 
42.33 
41.81 
41.41 
41.10 
40. 86 
40.67 
40.52 
40.40 
2.6421 2.8920 
3.4342 3.4361 
1.8288 2.1395 
3.3266 3.3457 
1.3339 1.6371 
3.0738 3.1061 
0.9763 1. 2692 
2.7515 2.7885 
0.7411 1.0033 
2.9132 2.4507 
0.5096 0.7541 
2.0859 221215 
0.3659 0.5797 
oalahiaay 1.8109 
0.2021 0.3995 
1.4970 1.5271 
0.1236 0.2882 
1.2422 1.2694 
0.0164 0.1688 
1.0174 1.0411 
-0.0197 0.0992 
0.8190 0.8399 
-0.0756 0.0302 
0.6483 0.6660 
-0.1571 -0.0082 
0.5022 0.5175 
-0.2468 -0.0406 
0. 3801 0.3926 


37, 


3.0750 
3.4366 


22-4055 
323519 


1.9113 
3.1168 


1.5381 
2.8009 


1.2473 
2.4632 


0.9836 
2.1334 


0.7809 
1.8220 


0.5862 
1.5371 


0.4447 
1.2784 


0.3114 
1.0491 


0.2142 
0.8469 


0.1318 
0.6719 


0.0706 
Qe227, 


0.0263 
0.3968 


e2.6.d 
Oale it 


2250.4 
orer .¢ 


Cpe. 
petit 
O6a2 5% 
7608 .. 


Htet 
S042,% 


aCe of 
rae 


eoere 
pice. 


-8t¢.U 
ifz#.! 


. 

Ta Thelen 
eT. 
oc at 
ehhh 
Loulr 
ert? 
TER 
oF .0? 
t/he 
Ls 4a 
Cee) 
(pO st® 
er ve 
cc cf 
‘ere 
teyte 
ht ute 
Ph ye 
Te ,ve 
ty ee 
oP ,oF 

Ustian’ | vena } 
DHF Oak Rete. 
ety. once 
tfe,F =v BME 
taal wets 
Paget erin,t 
a) (ote ft 
Anet.§ OTRT 
Erun.' 4. 
Hz7e ck Ate .S 
i ee>.0 aga? 
ear .a vee, .) 
reve. vral.3 
eore.t eer yd 
Perit bint 6 
TEFL etre .f 
yaar.4 adc .f 
PAL Ad Shec.r 
wes) .@ aero.d 
rae /f erro. 
see6,o TVhO.0- 
ePtw.0 Cerh.& 
PUL ea 
#0 ade o 


seetse” 


1av ee SS 


_ 


} 7 i 
SE 
Te hee 


a+ 

: 69° : 

pul tet ecansee eepsssuhe dhe ube 
. 


ad dead 
» * 
» o 
e = 
‘ i 
. S 
> © bd y 
> 7 : bal 
° sos 
* 1. os 
Le jf 
* . hd 
. - - 
» : * « 4 
* - ts 
4 . a1, 
* . a 


i 


wie ee P| 


x 
abe 


bete,s 
4G0est 


> 


res 


(iSey! 
ters.) 


4400.9 
serp.i 


14 To 4h 


_— os 


; 


Pri 
“0Fi. 


J 


leabak 


peer 0 
woul, ; 
vaett 
vote: 
bu Ag? 


a * ~ - 
Ste tk 


ei 


af90,0- 
rSwact 


a i 


a 


— 
— 


ot ete 


eed 
= 


7 _ 


eeranieyaiainiapaceamyasceat 


{Teseeurdl dons 
ee nore 


Sh m8 habbo aati 
miniaenen inn 


